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A  MODERN  CONTROL  DESIGN  METHODOLOGY 
WITH  APPLICATION  TO  THE  CH-47  HELICOPTER 

Richard  D.  Holdridge,  Ph.D. 

Stanford  University,  1985 

A  control  system  design  methodology  is  developed  which  produces  robust, 
low-order  optimal  controllers  for  multiple-input  multiple-output  systems.  The 
methodology  attempts  to  focus  the  strengths  of  recent  “Modern  Control"  design 
algorithms  on  the  problems  associated  with  real  control  system  designs.  The 
methodology  is  a  set  of  procedures  which  aids  the  engineer  in  creating  a  realizable 
controller  in  either  digital  or  analog  form. 

To  demonstrate  the  usefulness  of  the  methodology,  two  control  augmentation 
systems  (CAS)  were  designed  and  flight  tested  on  a  CH-47  helicopter  at  NASA 
Ames  Research  Center.  The  first  design  was  a  longitudinal  cruise  CAS  giving 
the  pilot  decoupled  control  of  forward  velocity  and  climb  rate.  This  design  task 
demonstrated  the  low-order  controller  and  robustness  features  of  the  methodology. 
It  also  demonstrated  the  use  of  modern  control  techniques  in  designing  integral- 
error  controllers.  Flight  test  results  are  presented.  The  second  controller  is  a 
translational  velocity  command/  precision  hover  hold  system.  This  two  mode  con¬ 
troller  demonstrates  the  methodology  as  applied  to  a  more  complicated  design  task 
which  includes  control  law  switching  and  inner  loop/  outer  loop  considerations. 
Flight  test  results  are  also  presented. 
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State  Equations 
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(2.2) 

Measurements 


'1  0  0  0  u  [0  o' 

0-10  176.0  w  +  0  0 

0  0  0  1  l  0  0 

L  J  L  9  J 


0  0 
0  0 


(2.3) 


'.318  0  0  ' 

R  =  0  .318  0 

0  0  .00039, 


(2.4) 

(2.5) 


where: 

u-  forward  velocity,  ft/sec 

w-  vertical  velocity,  ft/ sec 

q-  pitch  rate,  radians [ sec 

9-  pitch  angle,  radians 

elevator,  radians 
fit  -  throttle,  ft/ sec2 

longitudinal  gusts,  ft/ sec 
w„-  vertical  gusts,  ft/ sec 

Vy.vj'.vg-  measurement  noise,  ft/ sec.  ft/sec. radians 
disturbance  spectral  density  [2j 

measurement  noise  spectral  density( assumes  standard  deviations  of  1  ft/ sec  for 
u  and  h  and  2  deg  for  9  with  .16  second  correlation  times(T'c)).  The  spectral 
densities  come  from  the  approximation  S  D.  fts  2a2Tc.[9] 


Figure  2.1:  Navion  Linear  Model  Before  Scaling.  Nominal  airspeed  is  100  knots 
(100  ft/ sec).  Note  the  large  order  of  magnitude  differences  in  the  elements  of  the  system 
matrices.  This  is  caused  by  the  different  units  of  the  states. 
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needed  is  a  linear  system  in  state  variable  form: 

x  =  Fz  +  Gu  +  Tu> 
y  =  Hx  +  Lu  +  Nw  +  v 

(2.1) 

i  =  Az  4-  By 
u  =  Cz  +  Z)y 

where: 

z-  system  states,  n  x  1 
z-  compensator  states,  r  x  1 
u-  controls,  m  x  1 

plant  disturbances,  m'  x  1 
y-  sensor  measurements,  p  x  1 
v~  sensor  noise,  p  x  1 

Q-  plant  disturbance  spectral  density  matrix,  m'  x  m' 

R  sensor  noise  spectral  density  matrix,  pxp 

This  model  is  needed  for  all  plant  conditions  for  which  the  controller  is  expected 
to  operate.  Typically,  the  conditions  include  the  nominal  operating  point  and 
several  off-nominal  conditions  for  which  the  control  system  must  also  be  stable 
and  perform  adequately.  Figure  2.1  shows  the  Navion  model  at  the  sea  level,  104 
knot  flight  condition.  No  off-nominal  flight  conditions  are  included  in  the  example. 
With  the  dynamic  model  now  available,  the  methodology  can  continue  with  the 
scaling  of  the  model. 

2.2.  Model  Scaling 

Model  scaling  is  the  systematic  process  of  changing  the  units  of  the  variables 
in  the  plant  model.  This  process  is  described  and  used  by  Bryson(BR,sect  7.2) 
and  is  a  necessary  first  step  in  the  process  of  compensator  order  reduction.  Since 


Chapter  2. 

The  Design  Methodology 


From  the  discussion  of  Chapter  1,  the  need  for  a  structured  approach  to  using 
modern  control  methods  to  design  control  systems  is  apparent.  This  chapter  de¬ 
scribes  a  design  methodology  developed  to  take  advantage  of  the  strengths  of  the 
modern  control  techniques  while  eliminating  or  reducing  their  weaknesses.  The 
methodology  is  a  refinement  and  an  expansion  of  design  procedures  described  by 
Bryson  and  taught  in  the  advanced  flight  control  course  at  Stanford  University.fi j 
Figure  1 .3  shows  a  flow  chart  describing  the  methodology.  This  chapter  describes 
each  step  and  explains  why  it  is  necessary.  To  facilitate  the  description  of  the 
methodolgy,  a  simple  design  task  is  traced  through  the  process,  namely  the  design 
of  a  longitudinal  autopilot  for  commanding  aircraft  airspeed  and  climb  rate.  The 
plant  model  represents  the  Navion  general  aviation  aircraft  at  sea  level  and  about 
100  knots(176  ft/sec)  airspeed.  The  objective  is  a  dynamic  compensator  which 
could  be  implemented  on  a  digital  compvter  in  the  aircraft. 

2.1.  Model  Derivation 

The  development  of  dynamic  models  for  use  in  control  system  synthesis  is 
an  entire  engineering  discipline  in  itself.  For  this  design  methodology,  the  model 
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dures  is  the  lack  of  experience  in  real  world  application,  Chapters  3  and  4  present 
the  results  of  applying  these  techniques  to  a  real  world  problem.  The  “real  world” 
problem  is  the  design  of  a  control  system  for  a  Boeing- Vertol  CH-47  “Chinook” 
helicopter  operated  by  the  NASA  Ames  Research  Center.  Chapter  3  shows  the 
application  of  the  methodology  to  the  design  of  a  longitudinal  cruise  control  aug¬ 
mentation  system(CAS).  Specifically,  the  CAS  is  designed  to  give  the  pilot  a  veloc¬ 
ity/climb  rate  command  capability  for  cruise  flight.  The  controller  is  MIMO  with 
4  measurements  and  2  controls.  Chapter  4  applies  the  techniques  to  the  more  dif¬ 
ficult  task  of  designing  a  position  hold/velocity  command  controller  for  the  CH-47 
during  hover  and  low  speed  operations.  For  this  task  there  are  two  modes  of  oper- 
ation(velocity  command  and  position  hold),  switching  logic  for  going  between  the 
modes,  11  measurements,  and  3  controls  which  must  be  considered  in  the  design. 
For  both  these  designs,  simulation  and  flight  test  results  are  presented. 


Figure  1.3:  The  Design  Methodology.  The  many  loops  indicate  the  iterative 
nature  of  the  process.  The  computer  programs  associated  with  each  step  are  described  in 
the  appendices. 
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•  The  design  and  analysis  tools  for  LQG  techniques  have  only  recently  be¬ 
come  easy  to  use  in  practical  engineering  applications. 

These  disadvantages  provide  a  dilemma  to  the  engineer  faced  with  a  difficult 
MIMO  control  design  task.  If  he  chooses  to  use  LQG  methods,  he  may  be  break¬ 
ing  new  ground  when  he  finally  implements  his  system  in  actual  operational  hard¬ 
ware/software.  On  the  other  hand,  if  he  uses  classical  techniques  in  his  design 
procedure,  he  may  end  up  with  a  design  lacking  in  performance. 

This  report  presents  a  control  system  design  methodology  which  attempts 
to  minimize  the  disadvantages  of  many  modern  control  techniques,  as  described 
above,  while  taking  advantage  of  the  strengths  of  both  modern  and  classical  con¬ 
trol  techniques.  Specifically,  the  methodology  gives  the  engineer  a  structured 
approach  to  designing  arbitrary  order,  robust,  optimal  controllers  which  can  be 
implemented  in  digital  or  analog  hardware.  In  this  case,  “arbitrary  order”  implies 
feedback  compensation  of  any  size,  regardless  of  the  plant  model.  This  is  an  im¬ 
portant  consideration  for  high  order  plants.  “Robust”  is  defined  as  insensitivity 
of  the  controller  performance  to  changes  in  the  plant.  “Optimal”  relates  to  the 
design  methodology’s  attempt  to  minimize  a  quadratic  performance  index.  This 
methodology  is  a  refinement  and  expansion  of  techniques  described  by  Bryson  [2] 
and  depends  heavily  on  a  robust-control  design  algorithm  by  Ly.[3]  Figure  1.3  is 
a  flow  diagram  showing  the  steps  involved  in  the  methodolgy.  Chapter  2  of  this 
report  is  a  detailed  description  of  this  methodology.  To  clarify  the  description,  a 
simple  design  is  taken  through  all  steps  of  the  design  procedures  with  explanations 
for  design  decisions  made  during  the  process. 

Since  one  of  the  major  disadvantages  of  using  modern  control  design  proce- 
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Modern  (or  LQG)  design  techniques  are  characterised  by  simultaneous  closing 
of  all  loops  from  measurements  to  controls.  Figure  1.2  shows  the  form  of  this  type 
of  controller.  LQG  design  techniques  offer  several  advantages  to  the  designer  of 
MIMO  controllers: 

•  They  ensure  a  stable  design,  if  it  exists,  for  the  given  measurements  and 
controls. 

•  They  result  in  coordinated,  “graceful”  controllers.  In  particular,  the  con¬ 
trols  do  not  fight  each  other  to  satisfy  the  design  requirements. 

•  For  complicated  systems,  the  design  can  be  done  fairly  quickly. 

Although  these  advantages  make  modern  control  design  techniques  very  powerful, 

there  are  disadvantages  which  cannot  be  disregarded: 

/J  ‘ 

•  The  designs  are  often  not  robust  to  changes  and  uncertainties  in  the  plant 
model,  i.e.  they  are  to  “finely-tuned”  to  the  plant  model. 

•  For  complex  systems,  the  controller  designs  are  complicated,  making  im¬ 
plementation  difficult. 

•  Only  a  few  designers  have  developed  physical  intuition  for  selecting  perfor¬ 
mance  indices. 

•  There  is  little  physical  intuition  developed  concerning  the  loops  closed  by 
the  design  process. 

•  There  are  relatively  few  examples  of  operational  systems  designed  using 
these  techniques.  Practicing  design  engineers  are  therefore  hesitant  to  use 
these  techniques. 
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Figure  1.1:  Classical  Control.  Characterized  by  incremental  loop  closures. 
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Figure  1.2:  Modern  Control.  Characterized  by  simultaneous  loop  closures. 
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1.2.  A  Modern  Control  Design  Methodology 

Classical  control  techniques  are  routinely  applied  to  the  design  of  MIMO  sys¬ 
tems.  The  process  is  characterized  by  successive  loop  closing.  (Figure  1.1)  Some 
of  the  advantages  of  using  this  approach  are: 

•  The  designs  are  simple  with  clear  physical  intuition  for  each  loop  closed. 

•  For  uncomplicated  or  open  loop  stable  systems,  the  procedure  is  fairly  fast. 

•  There  is  a  large  experience  base  for  applying  these  methods. 

•  Numerous  design  and  analysis  tools  exist. 

•  Operational  hardware/software  implementation  is  well  understood  and  straight¬ 
forward. 

However,  there  are  also  some  disadvantages  to  using  classical  design  techniques 
for  MIMO  systems: 

•  For  complicated  or  highly  unstable  systems,  one  loop  may  not  stabilize  the 
system. 

•  It  may  be  difficult  to  decide  what  loops  should  be  closed.  This  is  espe¬ 
cially  true  in  complicated  systems(large  space  structures  for  example),  or 
in  unfamiliar  systems  where  little  physical  intuition  has  been  developed. 

•  Different  loops  can  “fight”  each  other  resulting  in  using  more  control  than 
is  necessary  or  available. 

•  Many  of  the  design  techniques(Root  locus,  frequency  response, etc.)  are 
difficult  to  apply  to  MIMO  systems. 
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•  Better  Numerical  Methods(HQR,  HQZ,  SVD,  etc.)-  1960-present 

•  Singular  Value  Techniques-  Doyle  and  Stein  and  others,  1970-present 

•  Gradient  Design  Methods-  Ly  and  others,  1980-present 

These  techniques  are  often  labeled  “modern  control*  where  the  previously  men¬ 
tioned  methods(Root  locus,  Frequencly  Response,  etc.)  are  often  called  “classical 
control”  techniques.  Perhaps  the  major  difference  between  the  two  is  that  classical 
techniques  are  best  suited  to  single  input  single  output(SISO)  dynamic  systems 
where  the  modern  control  techniques  are  equally  applicable  to  multiple  input  mul¬ 
tiple  output(MIMO)  systems. 

Virtually  all  control  systems  in  operational  use  today,  in  numerous  diverse  ap¬ 
plications,  have  been  designed  using  classical  control  techniques.  These  techniques 
have  been  preferred  over  the  modern  control  methods  because  most  of  the  appli¬ 
cations  are  fairly  simple  SISO  systems.  Even  if  the  systems  being  controlled  aren’t 
truly  SISO,  acceptable  SISO  approximations  can  often  be  made.  For  instance,  air¬ 
craft  motions  can  be  separated  into  longitudinal  and  lateral  motions  which  can  be 
approximated  by  fast  and  slow  linear  SISO  systems.  Design  methodologies  based 
on  classical  control  use  these  SISO  approximations  for  autopilot  designs.  As  we 
might  expect,  this  approach  results  in  less  performance  than  might  be  possible  had 
we  not  made  the  simplifying  assumptions  (necessary  to  use  the  design  techniques). 
This  procedure  is  completely  inappropriate  if  the  system  to  be  controlled  cannot 
be  adequately  approximated  by  several  uncoupled  SISO  systems. 
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engineers  had  the  opportunity  to  build  greatly  improved  control  systems.  How¬ 
ever,  the  existing  analytical  tools  were  difficult  to  apply.  More  systematic  anal¬ 
ysis  and  design  tools  were  needed  to  aid  in  developing  control  systems.  Such 
techniques  were  developed  in  the  early  part  of  this  century.  A  few  of  these  tech- 
niques( applicable  to  linear  systems)  are  listed  below: 

•  Nyquist  Stability  Analysis-  Nyquist,  1930*8 

•  Frequency  Response-  Bode,  1940’s 

•  Root  Locus-  Evans,  1940’s 

McRuer’s  Chapter  1  contains  an  excellent  summary  of  these  early  analytical  tech¬ 
niques  and  their  applications.[lj  All  these  techniques(plus  others)  provided  the 
theory,  and  the  signal  processing  devices  existed  to  design  sophisticated  control 
systems  which  were  installed  in  operational  systems  such  as  aircraft,  missiles, 
rockets,  ships,  etc.  The  requirements  of  the  military  during  World  War  n  pushed 
the  development  and  refinement  of  these  techniques  and  of  hardware  capable  of 
implementing  the  designs. 

During  the  1950’s,  the  tools  and  hardware  available  to  control  engineers  con¬ 
tinued  to  expand.  It  was  during  this  decade  that  the  digital  computer  became  a 
practical  tool  for  designing  and  implementing  control  systems.  With  the  digital 
computer,  the  control  design  engineer  had  much  more  signal  processing  capability. 
Analysis  and  design  techniques  soon  developed  which  took  advantage  of  the  digital 
computer.  A  few  of  these  techniques  are: 

•  Kalman  Filter-  Kalman,  1960 

•  Linear  Quadratic  Gaussian(LQG)  techniques-  Bryson  and  others, 1960-present 
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Introduction 

1.1.  Historical  Background  on  Control  System  De¬ 
sign 

Control  system  design,  as  an  engineering  discipline,  is  relatively  new.  It  has 
been  coupled  closely  with  the  development  of  measurement  devices(sensors),  elec¬ 
tronic  signal  processing  capability(analog  and  digital  hardware),  and  the  develop¬ 
ment  of  mathematical  techniques  for  the  design  and  analysis  of  dynamic  systems. 
Before  electronic  signal  processing,  ingenious  mechanical  devices  were  designed 
which  incorporated  simple  feedback  mechanisms.  “Simple”  here  implies  their  an¬ 
alytical  complexity;  sometimes  these  devices  could  be  quite  compicated  mechani¬ 
cally.  Some  simple  examples  of  these  mechanical  feedback  devices  are: 

•  Water  level  control  in  a  tank  using  float  valves 

•  James  Watts’  centrifugal  governor  for  steam  engines 

Conventional  techniques  for  analyzing  dynamic  systems  (Newtonian  dynamics, 
Lagrangian  dynamics,  etc.)  were  adequate  to  study  and  design  these  mechanical 
devices. 

With  the  development  of  electronic  signal  processing  and  electronic  sensors, 
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many  dynamic  systems,  including  the  example  here,  have  parameters  with  different 
units,  it  is  often  difficult  to  compare  control  gains  or  estimator  gains  for  states 
or  measurements  with  differing  units.  In  the  Navion  example,  to  determine  the 
relative  importance  of  velocity  to  elevator  and  pitch  angle  to  elevator  gains  would 
be  difficult  since  velocity  has  units  of  ft/ sec  and  pitch  angle  has  units  of  radians. 
Model  scaling  allows  the  engineer  to  select  a  new  set  of  units  which  result  in  the 
plant  model  having  variables  with  units  having  a  sort  of  “equivalent"  importance 
to  the  designer.  Of  course,  selection  of  such  a  set  of  units  requires  an  understanding 
of  the  dynamics  of  the  process  being  controlled.  For  the  Navion,  the  angle  variables 
are  changed  from  radians  to  .01  radians  while  the  linear  variables  remain  in  units 
of  ft/ sec.  Such  a  selection  is  reasonable  since  .01  radians  and  1  ft/ sec  are  similar 
in  importance  to  a  pilot  flying  the  aircraft.  For  a  supersonic  aircraft,  a  suitable 
selection  might  be  .01  radians  and  10  ft/ sec. 

Figure  2.2  shows  the  scaled  dynamic  model.  This  change  of  scaling  is  accom¬ 
plished  analytically  using  similarity  transformations  based  on  the  changes  in  units. 
Appendix  A  presents  the  derivation  of  the  transformations  and  the  equations  used 
to  transform  all  matrices  of  the  dynamic  model  into  the  new  units.  The  computer 
program  ROPTSYS,  described  in  Appendix  J,  implements  these  equations  as  the 
first  step  in  calculating  an  optimal  full  order  compensator.  The  next  step  in  the 
design  methodology  uses  this  scaled  dynamic  system  as  a  starting  point. 

2.3.  LQG  Design 

During  this  step  in  the  design  procedure,  the  optimal  regulator  and  estimator 
gains  are  calculated  for  the  scaled  dynamic  system.  Appendix  B  shows  the  deriva- 
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where: 

it-  forward  velocity,  ft/sec 

w-  vertical  velocity,  ft/ sec 

q-  pitch  rate,  .01  radians /sec 

0-  pitch  angle,  .01  radians 

St~  elevator,  .01  radians 

St-  throttle,  ft/sec 2 

«»-  longitudinal  gusts,  ft/sec 

w„-  vertical  gusts,  ft/sec 

Vti,v-h,vg-  measurement  noise,  ft/sec,  ft/sec,  .01  radians 
Q-  disturbance  spectral  density,  (From  Bryson  [2]  sect  9.4) 

R-  measurement  noise  spectral  density  (assumes  standard  deviations  of  1  ft/sec 

for  u  and  h  and  2  deg  for  9  with  .16  second  correlation  times(Te).  The  spectral 
densities  come  from  the  approximation,  S.D.  «  2<r2Tt-[9] 


Figure  2.2:  Navion  Linear  Model  After  Scaling.  The  results  of  the  scaling  are 
most  noticeable  in  the  F ,  G,  and  R  matrices. 
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tion  of  the  technique  used  to  calculate  these  gains.  The  optimal  gain  calculation 
of  this  appendix  is  based  on  the  eigenvector  decomposition  of  the  Hamiltonian 
matrix  as  described  in  Hall  and  Bryson.[17]  A  modification  to  the  technique  is  the 
capability  to  weight  the  linear  combination  of  states  and  controls,  y0  =  H0x  +  Du, 
in  determining  regulator  gains,  and,  by  duality,  to  include  plant  noise  in  the  mea¬ 
surement,  ym  =  Hmx  +  Lu  +  N w  +  v,  in  determining  estimator  gains.1  This 
modification  is  essential  for  aerospace  applications  since  it  allows  acceleration  to 
be  used  as  a  measurement  for  estimator  design,  and  weighted  in  regulator  design. 
The  ROPTSYS  computer  program,  described  in  Appendix  J,  calculates  these  reg¬ 
ulator  and  estimator  gains.  It  also  calculates  the  compensator  dynamic  system 
based  on  these  gains.  One  item  of  note  is  that  the  A  and  B  matrices,  used  in 
the  quadratic  performance  index(P./.  =  Aya  4-  uT Bu)  dt ),  can  start  as  the 

identity  matrices  since  the  system  has  been  scaled  in  “equivalent”  engineering 
units. 

The  compensator  dynamic  system  is  displayed  in  three  formats: 

Conventional  Form 


z  =  {F-GC-KH)z  +  Ky 
u  =  Cz 


(2.9) 


where: 

z-  compensator  states,  n  x  1 
y-  sensor  measurements,  p  x  1 
u-  controls,  m  x  1 

K-  steady  state  Kalman  filter  gains,  n  x  p 
C-  optimal  regualtor  gains,  m  x  n 

'Franklin  and  Powell  show  this  derivation  for  the  direct  digital  design  case. (6l 
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Modal  Form 


Z  =  F mod  2  +  Kmod  V 


V  =  Cmoi  Z 


(2.10) 


where: 

Fmo<r  modal  form  of  (F  -  GC  -  KH) 
Kmof-  modal  form  of  K 
Cmod~  modal  form  of  C 
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Block  Minimal  Form 


Z  —  F 7i  in  Z  ■(■  F m  in  J I 
^  =  Cmm  2 


(2.12) 


where: 

Fmt„-  minimal  form  of  (F  —  GC  —  KH) 
Km  in  ~  minimal  form  of  K 
Cm minimal  form  of  C 
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The  modal  form  is  calculated  using  the  eigenvector  matrix  as  a  similarity  transfor¬ 
mation.  Appendix  C  shows  the  equations  needed  to  do  this.  The  transformation 
to  minimal  form  is  derived  in  Appendix  C.  The  modal  and  minimal  realizations 
are  important  since  they  are  unique  realizations  of  the  dynamic  system.  From 
the  block  minimal  realization  of  the  compensator,  the  designer  can  assess  the  rel¬ 
ative  importance  of  the  different  modes.  To  aid  in  this  determination  of  modal 
importance,  the  ROPTSYS  program  also  calculates  two  input/output  measures, 
both  based  on  work  done  by  Bernard. [7]  One  measure  is  the  singular  value  of  the 
residue  matrix  associated  with  that  mode.  It  is  defined  as: 


Ml  =  o(Ri)  =  vW  +  gzgjsjhjhi  +  h^h2  (2.14) 

where: 

the  two  rows  of  /fm,„  associated  with  the  mode  being  analyzed 
hi,  h 2-  the  two  columns  of  Cm,„  associated  with  the  mode  being  analyzed 
ff(iZi)-  the  singular  value  of  the  residue  matrix  for  the  i,r>  mode 

The  second  measure  of  merit,  M2,  is  just  the  first  weighted  by  the  real  part  of  the 
mode  being  evaluated: 

7f(  RA 


[2.15 
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The  larger  these  measures  are,  the  more  important  that  mode  is  to  the  overall 
input/output  characteristics  of  the  compensator.  In  the  case  of  unstable  compen¬ 
sator  modes,  these  should  not  be  eliminated  since  they  are  usually  required  for 
stabilization  of  unstable  plant  modes.  Figure  2.3  shows  the  results  for  the  Navion 
control,  which  came  from  ROPTSYS.  The  input  data  required  to  get  these  results 
are  shown  in  Appendix  J. 

From  these  results,  we  note  that  the  open  loop  Navion  has  a  well  damped  short 
period  (the  fast  open  loop  mode)  but  the  phugoid  (the  slow  open  loop  mode)  is 
only  lightly  damped.  The  regulator  design  improves  these  dynamics  by  forcing 
the  two  phugoid  roots  onto  the  real  axis  and  speeding  them  up.  The  resulting 
compensator  is  fairly  fast  (as  fast  as  the  aircraft  response  itself)  and  well  damped. 
This  design  was  also  analyzed  using  the  RSANDY  program,  described  in  Appendix 
K,  to  evaluate  the  compensator  performance  in  the  presence  of  plant  disturbances 
(vertical  and  longitudinal  wind  gusts)  and  measurement  noise.  The  wind  gust 
data  come  from  a  turbulence  model  described  by  Bryson. [1]  The  performance  is 
tabulated  in  Figure  2.4. 

Using  these  data,  Figures  2.3  and  2.4,  we  can  proceed  to  the  next  step  of  the 
design,  simplification  of  the  full-order  compensator. 

2.4.  Compensator  Order  Reduction  and  Reopti¬ 
mization 

As  with  each  step  in  this  methodology,  the  subject  of  compensator  order  reduc¬ 
tion  and  simplification  is  an  engineering  discipline  in  itself.  For  this  methodology, 
the  compensator  mode  measures,  Mi  and  M2,  are  the  criteria  for  deciding  which 
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Figure  2.3:  Navion  Full  Order  Compensator.  The  unsealed  matrices  are  in  units 
of  the  physical  system. 
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Figure  2.4:  Navion  Full  Order  Compensator  Performance.  This  statistical 
performance  is  based  on  the  noise  characteristics  described  in  Figure  2.1. 
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modes  of  a  compensator  can  be  eliminated.  The  Af2  term  seems  to  be  most  useful 
in  the  work  here  on  flight  vehicles.  A  mode  can  probably  be  eliminated  with  little 
performance  impact  if  its  value  of  A/2/ is  .1  or  less.  When  this  ratio  is 
greater  than  .1,  that  mode  may  still  be  neglectable  but  the  decision  must  rest  on 
reduced  order  analysis  results.  In  addition  to  eliminating  modes,  a  compensator 
can  also  be  simplified  by  eliminating  unimportant  measurements.  This  is  desir¬ 
able  since  it  would  lower  the  cost  of  the  control  system.  The  scaled  ifm,„  matrix 
provides  the  data  which  are  used  to  decide  if  a  specific  measurement  is  neces¬ 
sary.  If  the  magnitudes  of  the  elements  of  a  specific  column  of  Kmin  are  smaller 
than  the  other  column  elements,  the  measurement  associated  with  that  column 
may  be  neglectable.  Of  course,  the  final  decision  must  be  based  on  the  simplified 
compensator  analysis  results. 

Applying  these  design  guidelines  to  the  Navion  example,  Figure  2.3,  we  note 
the  two  real  modes  have  an  M2  smaller  than  the  single  complex  mode  so  we 
eliminate  these  modes  and  the  associated  rows  of  Kmi„  and  columns  of  Cmm. 
Examining  the  values  of  scaled  Kmi„,  we  note  the  column  associated  with  the  9 
(pitch  angle)  measurement  is  smaller  than  the  other  so  we  try  eliminating  this 
measurement.  The  importance  of  the  scaling  process  is  evident  here  since  the  9 
column  elements  in  the  unsealed  Kmin  matrix  have  large  magnitudes  due  to  the 
radian  units.  Thus,  without  scaling,  the  9  measurement  would  have  seemed  more 
important  than  it  was.  The  compensator  is  now  second  order,  2  input,  and  2 
output,  a  significant  simplification  from  the  fourth  order,  3  input,  and  2  output 
full  order  compensator. 

The  next  step  is  to  analyze  and  reoptimize  the  reduced  order  compensator 
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Simplified  Compensator  before  Reoptimisation 
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Simplified  Compensator  after  Reoptimisation 
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Figure  2.5:  Navion  Reduced  Order  Compensator  Design  Results.  The 

reduced  order  compensator  was  unstable  before  the  reoptimization.  The  performance  is 
based  on  the  disturbance  and  noise  properties  of  Figure  2.1 


using  the  RSANDY  computer  program.  The  RSANDY  program,  described  in  Ap¬ 
pendix  K,  is  a  modified  version  of  a  design  and  analysis  code  written  by  Ly.[8]  The 
program  designs  robust,  low  order  compensators  using  a  gradient  search  technique 
based  on  a  quadratic  performance  index.  The  results  of  using  this  code  on  the 
Navion  simplified  compensator  are  shown  in  Figure  2.5. 

The  reduced  order  compensator  resulted  in  an  unstable  closed  loop  system. 
This  points  to  the  necessity  of  the  reoptimization  step  after  compensator  sim¬ 
plification.  The  utility  of  Ly’s  code  is  also  evident  here  since  it  allows  unstable 
initial  guesses  when  reoptimizing  the  compensator.  Another  important  feature, 
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not  used  in  this  example,  is  the  capability  of  including  several  plant  conditions  in 
the  optimization  to  ensure  compensator  stability  robustness  to  changes  in  plant 
parameters.  The  selection  of  weighting  matrices  in  the  quadratic  performance  in¬ 
dex,  used  by  the  RSANDY  program,  was  based  on  “Bryson’s  Rule*  as  applied  to 
the  scaling  parameters.  Briefly,  “Bryson’s  Rule*  suggests  that  the  outputs  and 
controls  be  weighted  by  the  inverse  square  of  their  scale  factors.  For  this  example, 
the  outputs  were  u(velocity)  and  A(climb  rate),  and  the  controls  were  $e(elevator) 
and  <5, (throttle),  with  units  of  ft/sec,  ft/sec,  radians ,  and  ft/sec2,  respectively. 
The  scaled  units  were  ft/sec,  ft/sec,  .01  radians,  and  ft/sec2.  The  performance 
index  is  then: 


where: 


P.I.  =  f  '  (yTQy  +  uTRu)dt 
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(2.20) 


For  more  complicated  problems,  these  weighting  terms  may  need  to  be  adjusted 
to  get  the  desired  performance.  The  input  file  for  RSANDY  used  to  get  the  results 
of  Figure  2.5  is  listed  in  Appendix  K. 

The  results  of  Figure  2.5  show  what  might  be  expected.  The  simplified  com¬ 
pensator  has  slightly  degraded  performance  in  that  both  the  aircraft  motion  and 
control  activity  is  slightly  larger  than  the  full  order  compensator.  These  statistical 
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performance  parameters  are  important  since  they  indicate  the  disturbance  rejec¬ 
tion  properties  of  the  design  but  we  recall  the  design  objective  was  a  velocity  and 
climb  rate  command  system  for  the  Navion.  Since  the  designs  thus  far  are  only 
regulators,  we  still  need  to  dv  sign  the  command  capability  for  the  controller. 

2.5.  Command  Inputs 

Generally,  control  systems  are  expected  to  do  more  than  regulate  the  outputs 
of  a  process  to  some  nominal  value.  We  need  the  ability  to  change  the  output 
of  the  compensated  system  to  any  selected  realizable  value.  To  accomplish  this 
end,  the  desired  outputs  are  used  to  determine  actuator  commands  which  are 
the  steady  state  controls  for  the  desired  outputs.  Figure  2.6  shows  this  concept. 
The  maximum  number  of  outputs  which  can  be  commanded  is  the  same  as  the 
number  of  independent  controls  in  the  system.  This  means  that  the  feedforward 
matrix  is  square.  An  additional  constraint  is  that  the  desired  outputs  must  be 
physically  realizable.  For  instance,  velocity  and  position  cannot  be  decoupled,  i.e. 
you  cannot  command  a  steady  velocity  while  holding  position.  Appendix  D  derives 
the  algorithm  which  calculates  this  decoupling  feedforward  matrix  for  the  case  of 
equal  numbers  of  controls  and  desired  outputs.  This  algorithm  is  implemented  in 
the  SETPNT  computer  program  described  in  Appendix  L. 

To  complete  the  Navion  autopilot  design,  we  need  to  include  the  climb  rate 
and  forward  velocity  command  capabilities  required  by  the  specifications.  For  this 
aircraft,  these  two  outputs  can  be  decoupled  (controlled  independently)  with  the 
two  controls,  elevator  and  throttle.  Using  the  SETPNT  program  with  the  Navion 
data(shown  in  Appendix  L),  the  feedforward  matrix  for  the  full  order  compensator 
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y«  -  Desired  Output 
Vm  —  Measurements  =[^^3 
uc  —  Command  Control-giSm/] 
ur  —  Regulator  Control  =  |?KthO 

ut  -  Total  Control-gKSig 
K  —  Control  Sensitivities 
N  —  Feedforward  Gains 


Figure  2.6:  Commanding  a  Desired  Output.  The  K  matrix  is  diagonal  and  based 
on  human  factor  considerations.  The  N  matrix  comes  from  an  analysis  of  the  system  in 
steady  state;  the  derivation  is  described  in  Appendix  D. 
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13: 


.002165  -.001578' 

.4505  .1398  . 


For  the  reduced  order  controller,  the  result  is: 


(2.21) 


,002243 

.5354 


-.002124 

.2653 


(2.22) 


where: 

<5e(rad»ans) 

.  6,(ft/sec 2) 

One  item  to  note  here  is  the  fact  that  the  feedforward  matrices  are  dependent  on 
both  the  plant  and  the  controller.  The  correctness  of  these  decoupling  feedforward 
matrices  becomes  clear  when  the  Navion  designs  are  evaluated  in  simulation. 


=  N 


Udt.irtd{ft/scc) 

hde.irtd{ft/sec) 


(2.23) 


2.6.  Simulation  and  Test 

The  final  step  of  any  design  methodology  is  the  demonstration  of  performance. 
The  first  step  in  such  a  demonstration  is  the  use  of  a  simple  linear  simulation.  Once 
the  design  is  validated  in  linear  simulation,  the  simulation  might  be  expanded  to 
model  some  of  the  important  nonlinearities  in  the  physical  system.  In  this  method¬ 
ology,  the  RSANDY  program,  described  in  Appendix  K,  includes  the  option  of 
creating  a  linear  simulation  model.  This  model  is  then  used  by  the  SIMPLOT 
program  (described  in  Appendix  M)  to  simulate  the  closed  or  open  loop  system 
with  or  without  the  random  disturbances.  Using  these  RSANDY  and  SIMPLOT 
programs  the  Navion  full  and  reduced  order  compensators  can  be  compared. 


The  results  of  step  commands  in  velocity  and  climb  rate  are  shown  in  the 
following  figures.  Comparing  Figures  2.7  and  2.8  we  note  the  performance  is 


9j 


Figure  2.7:  Velocity  Response  for  Full  Order  Compensator.  The  system  has 
a  rise  time  of  about  2.5  seconds  and  near  critical  damping. 


Figure  2.8:  Velocity  Response  for  Reduced  Order  Compensator.  The  rise 
time  is  about  2.5  seconds  but  the  damping  is  only  about  .5. 
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unsealed  model  q  and  9  appear  much  less  noisy  than  the  u  and  h.  measurements 
indicating  the  possibility  that  u  and  h  may  not  even  be  necessary.  After  scaling, 
the  spectral  densities  are  closer  in  magnitude  but  9  is  still  more  accurate.  With 
the  scaling  accomplished,  the  methodology  continued  with  the  design  of  the  full 
order  compensator  using  the  ROPTSYS  program. 

3.4.  LQG  Design  of  the  Full  Order  Compensator 


The  design  of  the  full  order  compensator  is  straightforward  since  we  have  all 
the  data  necessary  to  run  the  design  program  ROPTSYS.  Since  we  have  the  system 
scaled,  the  first  choices  for  output  and  control  weighting  matrices  in  the  regulator 
design  are  identity  matrices.  With  this  selection,  the  compensator  is  calculated  and 
shown  in  Figure  3.6.  Examining  the  closed  loop  roots  (the  estimator  and  regulator 
roots)  we  note  the  very  slow  estimator  pole  at  -.0088.  The  slow  estimator  pole 
results  from  having  a  very  accurate  measurement  of  9  (low  noise  spectral  density) 
but  relatively  noisy  measurements  of  u  and  h.  The  fast  filter  roots  estimate  mostly 
9  and  q  while  the  slow  mode  estimates  mainly  u.  If  we  had  introduced  h  and  u 
commands  into  the  compensator(Equation  3.9),  the  slow  estimator  mode  would 
not  have  affected  the  response  to  commands,  since  the  estimator  would  have  been 
“tracking”  closely  before  the  commands.  For  perfect  tracking,  the  estimator  modes 
are  not  excited  at  all.  This  need  to  include  the  command  in  the  compensator  was 
not  recognized  until  later  so  the  responses  of  Figure  3.7  show  the  poor  performance 
when  commands  are  not  properly  fed  forward  to  the  compensator.  Figure  3.8  shows 
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where: 

ix  forward  velocity,  ft/ sec 

w  vertical  velocity,  ft/sec 

7  pitch  rate,  .01  radians / see 

0  pitch  angle,  .01  radians 

ft  pitch  control,  .1  inches 

ft  collective  lever,  .1  inches 

tx»  longitudinal  gusts,  ft/ sec 

vertical  gusts,  ftjsec 

Vq.vs.v-f'.vu  measurement  noise,  .01  radians/ sec.  .01  radians,  ft/ sec.  ft/ sec 
Q  disturbance  spectral  density,  [l] 

R  measurement  noise  spectral  density(assumes  standard  devations  shown  in  Fig¬ 

ure  3.1  with  3  Hz  bandwidth)  The  spectral  densities  come  from  the  approxi¬ 
mation  S.D.  r*  2a2Tc  [9] 


Figure  3.5:  CH-47  Scaled  Longitudinal  Linear  Model  at  60  knots.  Even 
after  scaling,  the  9  measurement  is  much  more  accurate  than  the  others. 
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where: 

it  forward  velocity,  ft/ sec 

w  vertical  velocity,  ft/ see 

q  pitch  rate,  radians/ sec 

0  pitch  single,  radians 

fie  pitch  control,  inches 

hc  collective  lever,  inches 

u*  longitudinal  gusts,  ft/ sec 

w v  vertical  gusts,  ft/sec 

Vq,vg,vji,vu  measurement  noise,  radians/ sec.  radians,  ft/ sec.  ftj sec 
Q  disturbance  spectral  density,  [l] 

R  measurement  noise  spectral  density  (assumes  standard  devations  shown  in  Fig¬ 

ure  3.1  with  3  Hz  bandwidth)  The  spectral  densities  come  from  the  approxi¬ 
mation  S.D.  2a2  Te.  [9] 

Figure  3.4:  CH-47  Longitudinal  Linear  Model  at  60  knots.  The  highly 
accurate  9  measurement  comes  from  the  inertial  navigation  system(INS),  the  pitch  rate 
from  a  body  mounted  rate  gyro,  the  velocity  from  the  pitot-static  system,  and  the  climb 
rate  from  an  instantaneous  vertical  speed  indicator  (IVSI). 
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other  consideration  in  the  development  of  a  linear  model  for  this  design  is  the  need 
to  avoid  “nuisance  disengages”  of  the  experimental  control  system  in  the  CH-47. 

These  automatic  disengages  are  a  safety  feature  of  the  modified  control  system 
which  cause  the  experimental  control  system  to  be  tripped  off  when  the  control 
rates  exceed  a  known  rate.  The  approach  used  to  eliminate  the  possibly  of  these 
problems  was  to  include  simple  first  order  actuator  models: 

^actual  —  ^0  Ugdugj  +  40  Udciired  (3.2) 

By  weighting  uaetaai  in  the  performance  index,  the  control  rates  should  not  get 
large,  hence  eliminating  the  nuisance  disengages.1  Before  starting  the  control 
design,  the  linear  model  was  scaled  to  help  in  compensator  order  reduction. 

3.3.  Model  Scaling 

The  selection  of  units  for  the  CH-47  at  a  60  knot  flight  condition  is  dependent 

on  our  understanding  of  the  aircraft  dynamics.  In  this  case,  1  ft/ sec  of  velocity 

or  climb  rate  is  about  as  important  to  the  pilot  as  .01  radian  of  pitch  angle  or 

.01  radian/ sec  of  pitch  rate.  Similarly,  the  pilot  would  feel  comfortable  using  .1 

inch  of  either  control  to  command  the  1  ft /sec  of  velocity.  Another  way  of  looking 

at  this  scaling  is  to  note  that  a  .01  radian  pitch  angle  would  result  in  a  1  ft/ sec 

climb  rate  at  the  nominal  airspeed  of  60  knots(»  100/f/sec).  Figures  3.4  and  3.5 

show  the  linear  model  at  60  knots  both  before  and  after  scaling  the  system  into 

the  units  above.  1 

The  usefulness  of  the  scaling  step  is  again  evident  if  we  look  at  the  measure- 

ment  noise  spectral  density  matrix  for  the  scaled  and  unsealed  models.  In  the 

'Flight  test  later  confirmed  that  the  experimental  system  seldom  had  these  type  disengages.  | 
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yt  —  Actual  Outputs  (Subsat  of  Measurements) 
—  Desired  Outputs 
Xi  —  Measurements 
u«  —  Command  Control 
u!  —  Integral  Control 
ur  —  Regulator  Control 
—  Total  Control 


Figure  3.3:  Controller  Structure  for  The  Longitudinal  Cruise  Autopilot. 

Flexibility  comes  from  making  the  controls,  outputs,  and  measurements  vectors  rather 
than  scalars. 
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The  form  of  the  controller  to  be  used  for  the  design  depends  on  its  function. 
For  this  design,  the  controller  must  give  the  pilot  a  command  capability.  From 
Section  2.4,  we  recall  that  the  command  capability  comes  from  the  feedforward 
matrix  which  decouples  the  desired  outputs.  Appendix  D  shows  that  this  matrix 
results  from  inverting  the  steady-state  dynamic  system.  In  flight,  this  matrix  would 
exactly  decouple  the  two  desired  outputs  only  if  the  helicopter  were  actually  linear 
and  had  the  same  parameters  as  the  model.  Of  course,  this  can  never  happen  for 
several  reasons: 

•  The  helicopter  dynamics  are  not  linear. 

•  The  helicopter  has  dynamics  which  were  not  modeled. 

•  The  atmosphere  is  never  “standard”. 

•  The  helicopter  burns  off  fuel  during  flight,  changing  its  inertia  properties. 

•  The  sensors  have  biases  and  scale  factor  errors. 

The  list  could  go  on  but  it  is  clear  we  need  to  correct  for  differences  between  the 
model  and  the  actual  aircraft.  One  obvious  solution  is  to  use  integral-error  control. 
Figure  3.3  shows  a  form  of  the  controller  which  includes  this  capability. 

An  interesting  difference  between  this  type  of  controller  and  more  conventional 
flight  control  systems  is  the  presence  of  off-diagonal  terms  in  the  integral-error 
feedback.  Normally,  the  integral  feedback  is  from  desired  output  to  the  primary 
actuator  for  that  output.  For  instance,  a  classically  designed  integral  controller 
would  use  integral  of  velocity  error  for  pitch  control  and  integral  of  h  error  for 
collective.  These  off  diagonal  gains  would  be  difficult  to  find  using  classical  tech¬ 
niques  since  four  transfer  functions  would  be  used  to  calculate  the  four  integral 
gains.  Finding  these  gains  is  straightforward  using  the  RSANDY  program.  An- 
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Figure  3.2:  CH-47  Longitudinal  Model.  This  is  a  conventional  longitudinal  model 
found  in  most  textbooks.  The  climb  rate  measurement  is  based  on  the  approximation, 
k  —  to  H*  Unomi„al9. 

the  need  for  a  fairly  rigid  control  structure.  This  was  necessary  since  the  controller 
was  to  be  programmed  in  fixed  point  assembly  language  on  the  aircraft’s  Sperry 
18 19 A  flight  control  computer.  Since  major  controller  structural  changes  would  be 
difficult  and  time  consuming,  the  baseline  controller  had  to  have  a  broad  control 
structure  which  would  allow  design  flexibility  to  come  from  changes  in  the  com¬ 
pensator  order  or  in  the  matrices  themselves.  With  these  goals  and  constraints  in 
mind,  the  application  of  the  methodology  began. 


3.2.  Linear  Models  and  Basic  Control  Structure 

The  models  used  for  these  designs  come  from  Reference  11.  This  reference  gives 
the  aerodynamic  coefficients  at  several  different  flight  conditions  and  a  general  form 
of  the  “F”  and  “G”  matrices  based  on  these  coefficients,  the  inertia  properties  of 
the  aircraft,  and  the  nominal  airspeed.  Appendix  G  shows  this  general  form  as 
well  as  the  resulting  linear  models  used  for  this  research.  Figure  3.2  shows  the 
parameters  of  the  longitudinal  model. 
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Figure  3.1:  Sensors  Available  on  the  NASA  CH-47.  The  standard  deviations 
are  estimates  based  on  flight  test  data  from  the  CH-47  instrumentation  system.  The  .053 
second  correlation  times  corresponds  to  a  3  hz  bandwidth. 
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Similarly,  to  establish  a  constant  speed  climb,  the  pilot  must  add  collective  thrust 
to  begin  climbing  but  simultaneously  adjust  pitch  attitude  to  hold  airspeed.  In  a 
sense,  the  pilot  (through  his  extensive  training)  becomes  an  inner-loop  decoupling 
controller,  needed  to  give  good  speed  and  climb  rate  performance,  the  ultimate 
goal  in  cruising  flight.  The  design  goal,  in  terms  of  human  factors,  can  be  restated 
as  reducing  the  pilot’s  workload  in  cruise  flight  by  taking  him  out  of  the  aircraft 
inner  loops  and  giving  him  a  direct  velocity  command  system. 

With  the  design  goal  stated,  the  physical  hardware  and  flight  software  con¬ 
straints  must  also  be  considered.  Appendix  F  briefly  describes  the  highly  modi¬ 
fied  CH-47B  flight  research  helicopter  flown  and  maintained  at  the  NASA  Ames 
Research  Center.  The  Langley  Report  gives  a  more  detailed  description  of  this 
particular  helicopter,  which  has  been  modified  to  include  a  full  authority  fly-by¬ 
wire  flight  control  system  and  extensive  instrumentation.[10]  Figure  3.1  shows  the 
list  of  sensors  available  on  the  aircraft.  Unfortunately,  many  of  these  sensors  were 
quite  noisy  in  the  sense  that  they  contained  frequencies  associated  with  the  rotor 
motion.  Rather  than  attempt  to  use  the  noisy  measurements  directly  in  the  digital 
compensator,  the  TR-48  analog  computer  on  the  aircraft  was  programmed  to  filter 
several  of  the  important  aircraft  motion  sensor  outputs.  Specifically,  fourth  order 
Bessel  filters  were  used  on  the  longitudinal,  lateral,  and  vertical  accelerometers;  a 
third  order  Bessel  filter  was  used  on  the  roll  rate  gyro;  and  second  order  Bessel 
filters  on  pitch  rate,  yaw  rate,  altitude,  and  pitot-static  airspeed  measurements. 
The  filters  were  designed  with  5  Hertz  breakpoints  to  eliminate  the  “3  per  rev” 
and  “6  per  rev”  main  rotor  harmonics  at  11  and  22  Hz.  A  brief  description  of  these 
filters  is  included  in  Appendix  H.  A  further  constraint  on  the  design  process  was 
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Longtidinal  CAS  for  the  CH-47 

3.1.  Design  Goals  and  Constraints 

The  first  example  of  an  application  of  the  methodology  in  Chapter  2  is  the 
design  of  a  longitudinal  cruise  autopilot  for  the  CH-47.  The  goal  was  to  synthesize 
a  controller  which  gave  the  pilot  independent  control  of  airspeed  and  climb  rate 
using  separate  pilot  controls.  For  this  design,  the  pilot  longitudinal  stick  was  used 
for  the  airspeed  control  and  his  collective  lever  was  used  for  climb  rate  control. 
This  implementation  is  unconventional  since  most  helicopter  control  systems,  even 
the  highly  augmented  ones,  give  the  pilot  either  pitch  rate  or  pitch  angle  command 
from  longitudinal  stick  and  direct  collective  thrust  command  from  the  collective 
lever.  Some  of  the  more  modern  helicopters  (CH-47D,  HH60,  CH-53)  do  close 
outer  loops  such  as  altitude  or  speed-hold  around  the  collective  thrust  and  pitch 
rate  inner  loops,  but  these  are  usually  implemented  as  additional  pilot-selectable 
autopilot  modes.  One  trouble  with  pitch-angle,  collective-thrust  controllers  is  the 
coupling  between  the  two.  To  increase  speed  in  level  flight,  the  pilot  must  pitch 
down,  add  collective  thrust  to  maintain  altitude  and  accelerate,  then  pitch  up 
and  reduce  to  a  thrust  slightly  higher  than  the  original  to  hold  the  new  airspeed. 
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the  compensator  and  then  further  simplifies  the  discrete  compensator  by  putting 
it  back  into  block  minimal  form.  The  SETPNT  program  uses  an  algorithm  which 
finds  the  exact  digital  representation  of  a  2  x  2  compensator  subsystem.  One 
advantage  of  having  the  digital  compensator  in  block  minimal  form  is  the  reduction 
of  arithmetic  operations  required  by  the  computer.  Specifically,  the  minimal  form 
reduces  the  number  of  multiplies  and  adds,  required  by  the  compensator,  by  r(r—  1) 
where  r  is  the  order  of  the  compensator.  For  large  order  compensators,  such  as 
required  by  systems  controlling  structural  modes,  the  reduction  can  be  important. 
The  SETPNT  program  calculates  this  digital  compensator  in  a  form  immediately 
usable  in  a  floating  point  digital  controller. 

2.7.  Summary 

This  chapter  has  shown  an  approach  to  designing  control  systems  based  on 
“modern  control”  methodologies.  The  approach  is  especially  useful  for  MIMO  sys¬ 
tems  of  large  order  where  compensator  order  reduction  is  essential.  The  Navion 
autopilot  design  was  a  simple  example  chosen  to  illustrate  the  design  methodol¬ 
ogy.  In  applying  this  methodology  to  real  control  systems,  problems  will  surface 
requiring  the  designer  to  modify  and  iterate  the  design  to  achieve  the  desired  per¬ 
formance  specifications.  The  next  two  chapters  show  this  methodology  applied  to 
a  real  design  problem.  They  illustrate  the  usefulness  of  the  approach  applied  to  a 
fairly  complicated  design  task.  They  also  show  the  type  of  problems  which  surface 
in  real  design  applications. 
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Figure  2.12:  Navion  Reduced  Order  Compensator  Performance  in  Tur¬ 
bulence.  The  performance  is  nearly  identical  to  the  full  order  compensator  (Figure  2.11) 
and  is  shown  numerically  in  Figure  2.5 

could  be  “tuned”  by  changing  the  weighting  matrices  and  disturbance  properties 
in  the  RSANDY  program. 

Eventually,  the  control  system  may  be  put  into  the  physical  system  and  tested. 
Since  these  compensator  designs  can  be  computationally  intensive,  it  is  unlikely  to 
see  any  analog  or  continuous  implementations  for  anything  other  than  the  simplest 
designs.  Since  the  entire  design  process  has  been  in  the  continuous  domain,  the 
final  compensator  must  therefore  be  discretized  to  be  useful  in  a  digital  computer 
based  control  system.  The  process  of  discretization  is  simple  since  the  compensator 
is  in  a  block  minimal  form.  If  the  control  is  assumed  to  be  a  zero  order  hold 
(ZOH),  i.e.  step  commands  over  the  cycle  time,  then  the  discrete  form  can  be 
solved  exactly  for  the  2x2  and  lxl  blocks  which  make  up  the  compensator 
dynamic  system.  The  SETPNT  program  does  this  exact  ZOH  discretization  for 
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Figure  2.9:  Climb  Rate  Response  for  Full  Order  Compensator.  The  rise 
time  is  2  seconds  with  critical  damping. 


as  we  might  expect.  The  velocity  is  better  damped  and  faster  using  the  full 
order  compensator.  Also,  with  the  elimination  of  the  9  measurement,  we  note  the 
reduced  order  compensator  has  greater  pitch  angle  excursions  due  to  the  velocity 
command.  Comparing  the  climb  rate  responses  of  Figures  2.9  and  2.10,  we  note 
similar  performance  between  the  full  and  reduced  order  compensators.  To  compare 
the  disturbance  rejection  characteristics  of  the  two  designs,  we  use  the  SIMPLOT 
program  with  plant  disturbances  and  measurement  noise.  Figures  2.11  and  2.12 
show  the  system  performance  in  response  to  an  initial  pitch  rate  of  5  degrees  in 
the  presence  of  identical  disturbances.  As  expected  the  full  order  compensator 
does  a  better  job,  as  seen  earlier  in  Figure  2.5.  At  this  point,  the  design  could 
be  iterated  by  eliminating  more  modes  or  measurements,  or  by  adding  modes 
or  measurements  to  improve  performance.  Alternately,  the  reduced  order  design 
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the  responses  when  the  commands  are  correctly  accounted  for  in  the  compensator. 


z  —  4  z  ■+■  B  ym  4-  Tmjn  Tmoj  G  ueommanj 


M regulator 


=  Cz 


(3.9) 


where: 

A:  the  block  minimal  form  of  the  compensator  dynamics  matrix(F -GC-KH)  which 

already  includes  the  regulator  control  effects  (the  GC  term) 

B  :  The  minimal  form  of  the  Kalman  gain  matrix  from  ROPTSYS 

C  :  The  minimal  form  of  the  optimal  regulator  gain  matrix  from  ROPTSYS 

G:  the  control  distribution  matrix  for  the  physical  system 

r~‘n:  the  similarity  transformation  from  the  modal  form  to  the  block  minimal  form 
described  in  Appendix  C 

the  similarity  transformation  from  the  nominal  form  to  the  modal  form,  composed 
of  the  eigenvectors  of  F  -  GC  —  KH 

At  this  point  in  the  process,  alternative  approaches  to  the  design  were  used. (Still 
before  we  realized  how  to  implement  Equation  3.9  above.)  Two  such  alternatives 
are  shown  below: 

•  Pick  noise  spectral  density  matrices  independent  of  the  actual  noise  but 
which  give  good  time  responses. 

•  Use  an  inverse  optimal  solution  as  described  by  Bernard  to  establish  the 
performance  index  weighting  matrices  and  noise  spectral  density  matrices. [7] 

3.4.1.  Redesign  Using  Arbitrary  Measurement  Spectral  Den¬ 
sities 

The  first  intermediate  approach  to  improving  the  poor  time  response  in  the 
velocity  channel  required  selecting  noise  spectral  densities  which  give  an  adequate 
time  response.  This  approach  violates  the  assumptions  of  the  LQG  estimator 
design  procedure  where  the  Kalman  filter  gains  are  calculated  to  minimize  the 
estimate  errors  in  the  presence  of  the  given  noise.  However,  the  noise  properties  are 
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Velocity  Step  Climb  Rate  Step 


Figure  3.7:  Longitudinal  CAS  Time  Response!  using  Nomimal  Noise.  The 

poor  velocity  command  response  is  dominated  by  the  slow  estimator  mode  (s  =  -.0088). 


fairly  uncertain  and  furthermore,  the  goal  of  this  design  was  not  a  good  estimator, 
but  a  good  compensator.  The  selection  of  Q  and  R  matrices  which  give  “good” 
performance  is  an  iterative  process.  Figure  3.9  shows  the  compensator  when  the 
spectral  densities  of  all  the  measurements  are  set  to  .2.  The  very  slow  estimator 
pole  is  eliminated.  Figure  3.10  shows  the  improved  velocity  response  and  the 
unchanged  climb  rate  response. 

3.4.2.  Redesign  Using  an  Inverse  Optimal  Controller 

The  second  intermediate  approach  uses  Bernard’s  “inverse  optimal  controller” 
technique. [7]  This  technique  computes  AtB,Q  and  R  matrices  which,  if  used  in 
an  LQG  design,  give  a  closed  loop  system  with  desired  poles.  Figure  3.11  shows 
this  compensator  including  the  A,B,Q  and  R  provided  by  Bernard.  The  time 
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Figure  3.8:  Longitudinal  CAS  Time  Responses  with  Commands  to  the 
Compensator.  With  the  command  correctly  fed  to  the  compensator,  the  response  is 
identical  to  a  full-state  feedback  controller.  The  climb  rate  response  is  unchanged  from 
Figure  3.7. 
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Figure  3.9:  Longitudinal  CAS  using  Arbitrary  Measurement  Spectral 
Density.  The  slow  estimator  mode  of  Figure  3.6  is  eliminated. 
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Velocity  Step  Climb  Rate  Step 


Figure  3.10:  Longitudinal  CAS  Time  Responses  using  Arbitrary  Noise. 

The  poor  velocity  response  of  Figure  3.7  is  eliminated. 

responses  of  Figure  3.12  confirm  the  improved  performance. 

Figure  3.13  shows  the  statistical  performance  of  the  original  design  and  the 
two  alternate  designs.  As  expected,  the  nominal  design  has  the  lowest  errors.  The 
inverse  optimal  controller  is  nearly  as  good  but  the  other  is  clearly  the  worst. 
Based  on  the  combination  of  statistical  performance  and  the  time  responses,  the 
inverse  optimal  design  was  the  best  of  the  full  order  compensators.  With  the  full 
order  compensator  in  hand,  the  process  of  order  reduction  began. 

3.5.  Compensator  Order  Reduction 

Using  the  Af2  criterion  of  section  2.3,  both  the  nominal  and  the  inverse  optimal 
compensators  indicate  reduction  to  third  order  is  feasible.  From  a  more  practical 
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Mode 


Figure  3.11:  Longitudinal  Compensator  based  on  Inverse  Optimal  Solu¬ 
tion.  The  slow  estimator  mode  of  the  nominal  design,  Figure  3.6,  is  gone. 


ihhhmi 


Figure  3.13:  Performance  Comparison  of  Pull  Order  Compensators.  These 
data  are  based  on  the  measurement  noise  shown  in  Figure  3.1  and  disturbance  noise  of 
Figure  3.4. 
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standpoint,  using  output  feedback  of  the  measurements  is  also  a  possibility.  This 
is  appropriate  if  we  recall  from  section  3.1  that  all  the  measurements  were  being 
analog  filtered  before  entering  the  flight  control  computer.  The  next  two  sections 
show  these  two  designs. 

3.5.1.  Robust  Longitudinal  Third  Order  Compensator 

Using  the  nominal  compensator  reduced  to  third  order  as  a  starting  guess,  the 
RSANDY  program  was  used  to  optimize  the  low  order  compensator  using  its  ro¬ 
bust  design  feature.  The  following  flight  conditions  were  used  with  the  weightings 
shown. 


Condition 

Airspeed(knots) 

Climb  Rate(ft/min) 

Weighting 

Nominal 

0 

.8 

Off  Nominal  1 

60 

500 

.05 

Off  Nominal  2 

60 

-500 

.05 

Off  Nominal  3 

40 

0 

.05 

Off  Nominal  4 

80 

0 

.05  | 

Stability  is  ensured  by  the  RSANDY  program  for  each  of  these  flight  condi¬ 
tions.  Figure  3.14  shows  the  results  including  the  closed  loop  roots  at  the  nominal 
flight  condition.  There  is  still  a  slow  mode  (s  =  -.059)  which  is  shown  in  the  time 
responses  of  Figure  3.15  indicating  the  need  for  integral  control. 

When  the  integral  control  loop  is  added  to  the  third  order  controller,  the 
RSANDY  program  was  again  used  to  set  the  four  gains  in  the  K[  matrix  of  Figure 
3.3.  The  program  did  not  converge  with  only  the  four  gains  of  Kj  allowed  to 
vary  in  the  optimization.  By  releasing  both  the  Kj  and  the  compensator  gains, 
the  compensator  of  Figure  3.16  emerged.  Figure  3.17  shows  the  associated  time 
responses.  The  same  approach  was  used  to  set  the  integral  gains  for  the  full-order 
compensator  when  the  integral  control  loop  was  added.  Figures  3.18  and  3.19 
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Real  Part 
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Freq(rad/sec) 

Freq(Hz) 

-75.45 

0.0 

1.0 

75.45 

12.01 

-40.00 

0.0 

1.0 

40.00 

6.37 

-40.18 

0.0 

1.0 

40.18 

6.40 

-3.12 

0.0 

1.0 

3.12 

.50 

-.96 

1.36 

.58 

1.67 

.27 

-.96 

-1.36 

.58 

1.67 

.27 

-.059 

0.0 

1.0 

.059 

.0094 

-.42 

0.0 

1.0 

.42 

.068 

-.96 

0.0 

1.0 

.96 

.15 

Figure  3.14:  Longitudinal  Reduced  Order  Compensator.  Reducing  from  6'h 
to  3rrf  order  decreases  the  number  of  independent  gains  in  the  compensator  from  36  to  18. 
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Velocity  Step  Climb  Rate  Step 


Figure  3.15:  Longitudinal  Reduced  Order  Compensator  Time  Responses. 

The  velocity  response  is  still  quite  slow  but  climb  rate  is  adequate. 


show  the  resulting  compensator  and  time  responses. 

3.5.2.  Longitudinal  Output  Feedback  Controller 

Since  the  measurements  available  (qffuh)  are  nearly  the  same  2  as  the  states  of 
a  fourth  order  model,  the  gains  of  a  regulator  design  are  a  good  starting  point  in 
using  RSANDY  to  compute  output  feedback  gains.  In  fact,  using  these  gains  with 
no  further  optimization  gave  quite  impressive  performance  as  shown  in  Figures 
3.20  and  3.21. 

As  with  the  reduced  order  compensator,  an  integral  control  loop  is  necessary 

to  ensure  that  the  actual  commands  are  achieved  in  flight.  Using  the  Kf  gains 

from  the  reduced  order  design(Figure  3.16)  and  the  output  feedback  controller  of 
JTbe  fourth  order  model  has  states  (uwq$)  and  h  unam,nai$  -  w. 
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0.0  1.0  0.0 

Amtn  =  -229.5  -16.82  0.0 

0.0  0.0  -2.48 


-22.83  - 

-190.9 

-.07 

.059 

57.00 

3528. 

.044 

-.33 

-18.57  - 

-59.93 

-1.37 

-.67 

[0.0 

-.1 

-.1 

m,n  1.45 

-.071 

.041. 

‘  -.0054  .0098 

.-.00823  .047. 


Closed  Loop  Eigenvalues 


(3.28) 

(3.29) 

(3.30) 

(3.31) 


Figure  3.16:  Longitudinal  Reduced  Order  Compensator  with  Integral 
Control.  The  i4min,  and  Cm,„  matrices  are  not  the  same  as  in  Figure  3.14. 
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Velocity  Step  Climb  Rate  Step 


Figure  3.17:  Longitudinal  Reduced  Order  Integral  Controller  Time  Re¬ 
sponses.  Although  the  velocity  response  is  still  fairly  slow,  this  design  was  selected  for 
evaluation  in  flight. 
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Figure  3.18:  Longitudinal  Full  Order  Compensator  with  Integral  Control. 

A3  with  the  third  order  design,  the  A,  K,  and  C  matrices  are  different  from  the  full  order 
design  without  integral  control  (Figure  3.11). 
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Figure  3.31:  Third  Order  Compensator  Flight  Response  to  Climb  Rate 
Command.  Atmospheric  turbulence  masked  the  transient  portion  of  the  response  but 
steady  state  was  achieved  in  about  10  seconds,  similar  to  the  simulation  of  Figure  3.17. 
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Figure  3.30:  Pull  Order  Compensator  Flight  Response  to  Climb  Rate 
Command.  As  in  the  velocity  response  of  Figure  3.27,  the  h  command  response  was 
well  predicted  by  the  simulation  results  of  Figure  3.19.  This  response  also  shows  the 
decoupling  between  velocity  and  climb  rate. 
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Figure  3.29:  Output  Feedback  Compensator  Flight  Response  to  Velocity 
Command.  As  in  the  responses  of  Figures  3.27  and  3.28,  this  flight  response  agrees 
with  the  simulation  predictions.  In  this  case,  however,  the  overall  loop  gain  (final  gain  to 
the  actuator)  was  reduced  by  half. 
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Figure  3.28:  Third  Order  Flight  Response  to  Velocity  Command.  This 
reponse  to  a  command  of  approximately  10  /t/see  follows  closely  the  simulation  results 
of  Figure  3.17. 
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Figure  3.27:  Pull  Order  Compensator  Flight  Response  to  Velocity  Com¬ 
mand.  The  rapid  velocity  response  and  poorly  damped  pitch  angle  were  also  evident  in 
the  simulation  response  of  Figure  3.19. 
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in  a  realistic  task  such  as  flying  a  precision  approach. 

3.8.  Summary  of  Results  for  the  Longitudinal  CAS 
Design 

The  results  of  this  section  can  be  separated  into  two  categories: 

•  the  design  task 

•  the  flight  test  results 

The  design  task  was  important  because  it  established  an  experience  base  for  use 
of  the  design  methodology.  This  task  showed  that  use  of  integral-error  feedback 
in  a  state  variable  based  controller  eliminates  the  effects  of  an  inaccurate  plant 
model  in  achieving  commanded  outputs.  It  also  showed  the  difficulty  encountered 
when  using  LQG  design  techniques  with  specified  time  domain  properties.  In 
this  example,  three  methods  of  setting  weighting  matrices  and  spectral  density 
matrices  were  compared.  From  a  practical  standpoint,  this  first  design  task  was 
very  important  since  an  accurate  way  of  scaling  the  analytical  design  for  use  in 
the  fixed  point  flight  computer  was  developed.  An  equally  important  result  of  this 
task  was  the  coding  and  testing  of  a  general  form  of  a  modern  controller  for  the 
Sperry  1819  flight  computer. 

The  flight  test  results  showed  the  technical  feasibility  of  decoupled  control 
using  a  “modern”  controller.  The  test  system  had  adequate  handling  qualities 
but,  due  to  time  constraints,  no  attempt  was  made  to  iterate  the  designs  for 
“good”  handling  qualities.  Also,  the  disturbance  rejection  capabilities  were  not 


thoroughly  investigated. 


Section  3.7 


62 


took  place  at  the  Crows  Landing  test  facility(in  the  San  Joaquin  valley)  to  avoid  the 
heavy  traffic  of  the  South  San  Francisco  Bay  Area  near  Ames  Research  Center.  The 
purpose  of  the  flight  test  was  to  validate  the  performance  of  the  different  controllers 
rather  than  “tune”  the  system  for  maximum  pilot  acceptance.  Step  commands 
from  the  computer  or  the  pilot  were  used  to  evaluate  the  different  systems.  The 
nominal  airspeed  was  60  knots  but  stability  and  performance  were  checked  from 
40  to  80  knots.  Figures  3.27,  3.28  and  3.29  show  the  responses  of  the  full-order, 
third-order,  and  direct  feedback  compensators  to  velocity  step  commands  of  about 
10  ft/ sec.  All  three  of  these  systems  show  good  decoupling  between  velocity  and 
climb  rate.  Their  velocity  responses  are  similar  to  the  simulation  results  of  figures 
3.12,  3.17,  and  3.23  but  the  pitch  angle  behavior  is  less  damped  than  the  simulation 
results.  Chen  has  shown  that  unmodeled  dynamics,  especially  rotor  dynamics,  are 
the  probable  cause  for  the  lower  achievable  control  bandwidth  in  flight.[12]  In  fact, 
the  output  feedback  flight  implementation  had  the  overall  pitch  and  collective  loop 
gains  reduced  by  half  to  achieve  an  acceptable  response.  Figures  3.30,  3.31  and 
3.32  show  similar  results  for  climb  rate  commands.  The  importance  of  integral 
control  is  demonstrated  in  figure  3.33  which  shows  the  response  to  a  climb  rate 
command  for  the  third  order  controller  without  integral  control.  Although  a  steady 
climb  rate  is  achieved,  an  unwanted  velocity  change  of  similar  magnitude  is  also 
present. 

Pilot  opinion  was  mixed  concerning  the  system.  They  were  impressed  at  how 
well  it  held  airspeed  and  climb  rate  but  the  transient  behavior  was  not  totally 
acceptable  and  the  system  was  “sloppy”  in  response  to  gusts.  This  is  valid  criticism 
based  on  the  pitch  angle  responses.  Unfortunately  the  system  was  not  evaluated 
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Figure  3.26:  Comparison  of  Analytical  Design  and  Flight  Test  Implemen¬ 
tation  of  the  Controller.  This  figure  confirms  the  correctness  of  the  steps  needed  to 
go  from  the  analytical  design  to  flight.  These  responses  are  for  the  3rd  order  compensator 
with  integral  control. 


tinuous  designs  were  digitized  and  scaled  correctly.  The  OBS  was  also  useful  in 
preliminary  pilot  evaluation  of  the  control  laws  and  in  initial  setting  of  the  stick 
sensitivities  and  other  pilot-related  items.  Figure  3.26  shows  a  comparison  of  the 
time  responses  for  a  climb  rate  command  for  the  third  order  system.  This  con¬ 
firmed  the  accurate  digitization  and  scaling  of  the  analytical  design.  All  flight 
controllers  were  similarly  checked  prior  to  each  flight. 


3.7.  Flight  Test  Results 


Once  the  designs  had  been  implemented  and  checked  out  in  the  flight  hardware, 
the  flight  testing  began.  Since  the  controller  was  designed  for  cruise,  the  testing 
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Figure  3.25:  Flight  Test  Implementation.  Sensor  outputs  and  actuator  commands 
were  filtered  in  the  analog  computer  while  the  control  laws  were  executed  in  the  digitial 
computer. 

pose  matrix  multiply  routine  was  programmed  to  take  advantage  of  the 
minimal  form  of  the  compensator  dynamics  matrix. 

•  The  existing  instrumentation  output  subroutines  were  modified  to  send 
internal  compensator  data  to  the  ground  support  station. 

•  Since  the  Sperry  1819A  flight  computer  is  an  18  bit  fixed  point  machine,  the 
matrices  from  the  previous  section  were  scaled  to  avoid  numerical  overflow 
during  program  execution.  Appendix  E  describes  this  “fixed  point  scaling” 
technique  and  lists  the  SCALEM  computer  programs  which  accomplish 
this. 

The  actual  assembly  code  implementing  the  controller  is  shown  in  Appendix  I. 

One  important  capability  of  the  research  system  was  onboard  simulation(OBS). 
The  OBS  allowed  the  real-time  flight  software  to  be  checked  in  the  closed-loop 
system  prior  to  actual  flight.  It  was  especially  useful  in  confirming  that  the  con- 
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could  work  without  integral  assistance.  This  would  determine  how  well  the  design 
models  compared  with  the  real  aircraft.  Using  the  SETPNT  program,  the  digi¬ 
tal  block  minimal  forms  of  each  of  these  compensators  were  computed  for  flight 
implementation. 

3.6.  Flight  Test  Implementation 

Once  the  analytical  designs  were  complete,  the  tedious  task  of  actual  flight 
implementation  began.  Figure  3.25  shows  a  block  diagram  of  the  control  structure 
on  the  research  vehicle.  The  TR-48  analog  computer  was  used  to  filter  the  sensor 
outputs  and  the  digital  actuator  commands  from  the  Sperry  computer.  The  digital 
commands  from  the  computer  were  filtered  to  avoid  possible  actuator  wear  caused 
by  the  20  Hz  chatter.  The  control  laws  were  implemented  in  the  Sperry  1819A 
digital  flight  control  computer.  Before  the  designs  of  the  previous  sections  could 
be  tested,  the  flight  control  computer  software  had  to  be  modified  to  use  them. 

In  order  to  implement  design  changes  more  quickly,  the  flight  software  was  set  up 
to  use  the  compensator  matrices  directly.  Some  of  the  considerations  involved  in 
programming  the  general  form  of  the  compensator  in  assembly  language  (shown 
in  Figure  3.3)  are  listed  below: 

•  Since  the  design  was  based  on  a  linear  perturbation  model,  the  measure¬ 
ment  and  control  trim  values  were  approximated  as  the  values  at  engage 

time;  these  were  subtracted  from  their  sensed  values  for  use  by  the  con-  , 

troller.  This  had  the  added  advantage  of  eliminating  engage  transients. 

•  The  compensator  states  were  initialized  to  zero  before  system  engage. 

•  To  minimize  time  required  for  the  compensator  calculations,  a  special  pur-  • 
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Velocity  Step 


Climb  Rate  Step 


Figure  3.23:  Longitudinal  Output  and  Integral  Feedback  Controller  Time 
Responses.  This  design  was  selected  for  flight  evalualtion. 
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Figure  3.24:  Performance  Comparison  of  All  Compensators.  The  boldbaced 
systems  were  selected  for  flight  evaluation. 
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Figure  3.22:  Longitudinal  Output  and  Integral  Feedback  Controller.  As 

with  the  full  order  and  3rd  order  designs,  the  use  of  integral  control  reduces  damping 
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Velocity  Step  Climb  Rate  Step 


Figure  3.21:  Longitudinal  Output  Feedback  Controller  Time  Responses. 

These  responses  are  critically  damped  and  typical  of  full-state  feedback  designs. 


Figure  3.20,  good  results  were  achieved  with  no  further  optimization.  Figures  3.22 
and  3.23  show  these  results. 

3.5.3.  Summary  of  Design  Results 

A  number  of  controllers  have  been  presented  to  show  the  iterative  nature  of 
the  design  methodology  and  to  show  an  application  of  several  methods  to  meet¬ 
ing  the  requirements.  Selecting  designs  for  flight  test  implementation  required  a 
review  of  these  results.  Figure  3.24  shows  a  comparison  of  all  the  designs.  The 
controllers  shown  in  boldface  in  Figure  3.24  were  selected  as  candidates  for  flight. 
The  three  designs  with  integral-error  control  show  a  comparison  of  three  sizes 
of  compensators  (full-order,  reduced  order  and  output  feedback).  The  3rd  order 
controller  without  integral  control  was  included  to  see  if  the  decoupling  matrix 
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Velocity  Step  Climb  Rate  Step 


Figure  3.19:  Longitudinal  Pull  Order  Integral  Controller  Time  Responses. 

This  design  was  also  selected  for  flight  evaluation. 


-5.74  -12.51  .041  .103 

.-1.23  -10.44  -.057  .013. 


(3.36) 


_ System _ 

Open  Loop(F) 
Closed  Loop(F  -  G  D  H„ 


_ Eigenvalues _ 

-40.0,  -40.0,  -2.54,  .503,  -.105  ±  276 
-39.4,  -37.52,  -2.07  ±  ;'.213,  -.72,  -.47 


_ Damping _ 

1, 1, 1,  unstable,  .35 
1,1,.  99, 1,1 


Figure  3.20:  Longitudinal  Output  Feedback  Controller.  Since  the  measure¬ 
ments  (q,  9,  h ,  and  u)  are  almost  the  same  as  the  states  of  a  4th  order  aircraft  model  (u, 
w.  q,  and  theta),  the  output  gains  were  set  to  the  full-state  feedback  gains  from  a  4,h  order 
regulator  design. 
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Figure  3.32:  Output  Feedback  Compensator  Flight  Response  to  Climb 
Rate  Command.  A3  in  Figure  3.29,  the  overall  loop  gains  to  the  actuators  were 
decreased  by  half  from  the  analytical  design. 
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Figure  3.33:  Third  Order  Flight  Response  to  Climb  Rate  Command  with¬ 
out  Integral  Control.  This  figure  documents  the  effectiveness  of  the  integral  control 
loops.  It  also  emphasizes  the  need  for  better  models  of  the  CH-47. 


Chapter  4. 


Hover  Controller  for  the  CH-47 


4.1.  Design  Goals  and  Constraints 

With  the  experience  gained  in  the  design  of  the  longitudinal  CAS,  a  more 
difficult  task  was  selected  for  the  next  flight  experiment.  The  second  application 
of  the  methodology  was  the  design  of  a  translational  velocity  command,  precision 
hover  hold  control  system  for  the  CH-47.  This  system  was  to  provide  the  pilot 
with  “split-axis”  control  of  translational  velocity  in  an  heading-oriented  inertial 
coordinate  frame.  “Split-axis”  here  means  the  pilot  could  select  either  a  velocity 
command  or  position-hold  control  mode  in  each  of  the  three  translational  degrees 
of  freedom  of  the  aircraft.  As  in  Chapter  3,  this  design  required  decoupling  of  the 
three  axes  of  interest.  Figure  4.1  shows  the  coordinate  system  used  for  the  design. 
As  in  Chapter  3,  this  control  law  is  somewhat  unconventional  since  the  pilot’s 
workload  is  reduced  by  removing  him  from  inner  loop  attitude  control  tasks.  For 
this  system,  forward  velocity  is  controlled  by  longitudinal  stick  displacement,  side 
velocity  by  lateral  stick,  and  vertical  velocity  by  the  collective  lever.  If  a  particular 
pilot  control  is  within  a  small  distance  of  the  neutral  position  (the  detent),  then 
the  system  enters  a  position-hold  mode  for  that  axis.  Yaw  rate  control,  using  the 
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Figure  4.1:  Hover  Control  System  Coordinate  System.  The  x ,  y,  and  z 

coordinates  are  in  a  heading-oriented  inertial  system. 


standard  CH-47  SAS,  comes  from  the  pedals.  No  operational  helicopter  control 
systems  have  this  type  of  capability  although  the  concept  was  tried  on  a  modified 
CH-47  helicopter  during  early  work  on  the  Army’s  Heavy  Lift  Helicopter(HLH) 
concept. [13]  This  type  of  control  law  has  numerous  advantages  with  potential 
applications  such  as: 


•  Search  and  rescue 


•  Shipboard  operations 

•  Slung  load  operations 

Before  proceeding  with  the  description  of  how  the  design  methodology  was 
applied,  the  design  constraints  must  be  mentioned.  As  in  the  longitudinal  CAS 
design,  the  TR-48  was  used  to  filter  the  aircraft  motion  sensor  data.  The  inertial 
velocity  and  position  measurements  came  from  either  a  laser  or  radar  tracker  on 
the  ground.  This  was  necessary  since  the  onboard  inertial  navigation  system(INS) 
had  a  drift  of  several  knots  (»  5ft/ sec)  and  provided  only  position  data  (no 
velocities).  To  avoid  major  changes  in  the  flight  software,  a  general  form  of  the 
controller  was  programmed  in  the  flight  computer  which  provided  for  a  flexible 
control  structure. 
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Figure  4.2:  CH-47  Hover  Model.  This  is  the  conventional  8,fc  order  model  described 
in  Appendix  G. 


4.2.  Linear  Model  and  Basic  Control  Structure 


The  8'*  order  basic  airframe  model  is  described  in  Appendix  G.  This  model  was 
augmented  with  three  actuator  states  as  shown  in  Equation  3.2.  These  actuator 
states  allow  the  the  designer  to  penalize  control  rate,  a  necessary  step  to  avoid 
nuisance  disengages  of  the  research  control  system  during  flight.  Figure  4.2  shows 
this  model  as  it  was  used  for  the  control  law  synthesis.  One  item  to  note  is  that 
there  is  no  yaw  control  in  the  model.  The  yaw  SAS  was  approximated  by  adding 
-1  to  the  Nt,  element  of  the  “F”  matrix  of  Figure  4.2  and  Appendix  G. 
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Since  the  longitudinal  CAS  showed  the  necessity  of  integral-error  control,  the 
hover  control  structure  had  to  accommodate  this  capability.  Figure  4.3  shows 
the  hover  controller.  Unlike  the  longitudinal  CAS,  integral  control  decoupling 
is  not  done  in  the  inner  loop.  Instead,  PID  (Proportional  Integral  Derivative) 
outer  loops  are  closed  separately  to  xeommanjf  y  command*  <md  zeommanj ,  which  act  as 
controls  for  these  outer  loops.  The  velocity  command  system  is  the  inner  loop. 
Setting  up  the  control  structure  in  this  way  had  several  advantages: 

•  The  inner  loop  and  outer  loop  designs,  both  fairly  complicated,  could  be 
separated. 

•  By  using  the  desired  outputs  as  controls,  the  magnitude  of  the  outer  loop 
gains  became  physically  meaningful.  This  proved  to  be  important  later 
when  these  gains  were  adjusted  during  flight  to  achieve  good  performance. 

•  Keeping  the  inner  loop  separate  made  the  mode  switching  between  pilot 
velocity  command  and  PID  control  easier. 

Before  the  full-order  design  began,  the  system  was  scaled  in  equivalent  units. 

4.3.  Model  Scaling 

The  hover  model  was  scaled  into  the  same  units  as  used  in  the  longitudinal 
CAS  design.  Ft/ sec  remained  ft/sec ;  radians  and  radian/sec  became  .01  radians 
and  .01  radians /sec]  and  inches  of  control  became  .1  inches  of  control.  Figures 
4.4  and  4.5  show  the  model  before  and  after  scaling.  The  ROPTSYS  computer 
program  did  this  scaling. 
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State  Equations 
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-0.021 

-0.00085 

0.0328 

-0.00019 

-0.137 

0.00265 

0.0248 

0.00374 

-0.296 

-0.00013 

-0.00852 

0.00058 

0.00925 

0.00017 

0.00234 

0.00039 

-0.00112 

0.00027 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0205 

2  585 

-0.106 

-31.986 

-1.494 

C.00414 

-0.165 

0.0292 

0.0419 

0.435 

0.362 

-3.71 

-0.716 

0.0382 

-0.0708 

0.0 

0.0427 

-1.23 

-0.00433 

0.0 

-0.0544 

-0.158 

-1.0 

0.0 

0.0 

1.0 

0.00788 

0.0 

1.0 

0.00091 

0.116 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.114 

0.939 

0.0 

31.99 

0.0118 

0.0635 

1.159 

0.0 

0.303 

-8.062 

.00002 

0.0 

-0.00596 

-0.0142 

.432 

0.0 

0.329 

0.019 

0.0 

0.0 

0.0461 

-0.00037 

.0425 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-40.0 

0.0 

0.0 

0.0 

0.0 

-40.0 

0.0 

0.0 

0.0 

0.0 

-40.0 

o 

o 

0 

1 .293  X  10"‘ 
0 
0 
0 
0 


0 

0 

0 

0 

3.232  X  10"5 
0 
0 
0 


0 

0 

0 

0 

0 

8.079  X  10"6 
0 
0 


0 

0 

0 

0 

0 

0 

3.232  X  10'7 
0 


0 

0 

0 

0 

0 

0 

0 

3.232  X  10'7 


Figure  4.4:  CH*47  Hover  Linear  Model.  The  last  3  states  are  40  radian  actuator 
models. 
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Figure  4.5:  CH-47  Scaled  Hover  Linear  Model.  After  scaling,  the  weak  coupling 
terms  in  the  F  matrix  become  very  obvious. 


Section  4.4 


78 


4.4.  LQG  Design  and  Compensator  Order  Reduc¬ 
tion 

The  scaled  11"*  order  model  shown  in  Figure  4.5  was  used  by  the  ROPTSYS 
program  to  calculate  the  full-order  compensator.  At  this  point,  an  important 
simplification  should  be  emphasized.  The  system  was  to  control  x,  y,  and  z  which 
are  inertial  velocities.  The  model  as  used  included  u,  v,  and  w  which  are  body 
axis  airmass  velocities.  The  implicit  assumption,  needed  to  facilitate  the  design, 
was  that  the  two  sets  of  velocities  were  equal: 

x  =  u 

y  =  v  (4.2) 

z  ~  w 

This  assumption  is  reasonable  only  if  6  and  <f>  remain  small,  which  they  must  for 
safe  hover  in  a  large  helicopter  such  as  the  CH-47.  The  outputs,  y0,  weighted  in 
the  performance  index: 

P  I-  =  f  {yj Ay„  +  uT Bu)  dt  (4.3) 

were  the  three  velocities  (u,  v,  and  u;)  and  the  three  control  rates.  Since  the 
system  had  been  scaled,  the  weighting  matrices,  A  and  B,  were  just  6x6  and 
3x3  identity  matrices.  Using  these  assumptions  and  criterion,  the  resulting  full- 
order  compensator  is  shown  in  Figure  4.6. 

Based  on  these  data,  especially  the  modal  cost  M2,  the  11'*  order  compensator 
was  reduced  to  5'*  order.  This  represents  a  significant  reduction  of  complexity  in 
the  resulting  compensator.  The  11'*  order  design  had  121  independent  gains  while 
the  5'*  order  compensator  had  only  55.  After  the  order  reduction,  the  compensator 
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Figure  4.6:  Hover  Full  Order  Compensator  Design  Results.  The  open  loop 
has  two  unstable  modes  which  are  divergent  pitch  and  roll  oscillations.  Figure  4.7  is  the 
reduced  and  reoptimized  compensator  based  on  this  full  order  design.  Order  reduction 
was  based  on  the  terms  shown  here. 
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was  optimized  using  the  RSANDY  program.  Figure  4.7  shows  the  optimal  reduced 
order  compensator  including  the  closed  loop  roots.  This  optimization  step  was 
difficult  since  the  initial  guess  (the  original  lllfc  order  compensator  reduced  to  5‘fc) 
was  quite  unstable.  This  caused  numerical  overflows  on  the  VAX  computer  used 
to  run  the  RSANDY  program.  Convergence  to  a  stable  5th  order  compensator  was 
finally  achieved  after  numerous  iterations  of  the  outputs  and  the  output  weightings 
of  Equation  4.3.  Figure  4.7  shows  these  outputs,  y„,  and  the  elements  of  the 
diagonal  A  and  B  matrices.  This  figure  also  shows  another  aspect  of  the  difficulty 
of  this  optimization.  The  entire  C  matrix  was  allowed  to  vary  which  meant  that 
there  were  60  gains  being  adjusted  by  the  RSANDY  program,  5  more  than  the  55 
independent  gains  of  a  minimal  realizaton. 

The  simulation  step  responses  are  shown  in  Figures  4.9,  4.10,  and  4.11.  All 
three  velocity  responses  look  good  but  the  pitch  angle  damping  has  several  over¬ 
shoots.  A  modal  analysis  later  confirmed  this  by  showing  the  mode  at  -.50±j2.03 
of  Figure  4.7  to  be  strongest  in  9.  With  the  velocity  inner  loop  regulator  designed, 
the  feedforward  matrix  was  calculated  for  direct  command  of  x,  y,  and  z.  The 
PID  outer  loop  gains  could  now  be  designed. 

4.5.  Outer  Loop  Design 

As  discussed  in  Section  4.2,  the  velocity  command  inner  loop  and  the  PID 
outer  loop  were  separate.  With  the  inner  loop  velocity  controller  set,  the  outer 
loop  design  began.  This  approach  (inner  loop  first  then  outer  loop)  is  common  in 
the  development  of  operational  aircraft  autopilots  except  that  the  inner  loops  are 
normally  designed  classically  using  incremental  loop  closure.  Initially,  these  PID 
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!  Performance  Index  Data  ! 

Outputs  or 
Controls 

Units 

Weighting 

u 

ft/ tec 

5  x  10* 

V 

/t/iec 

5  x  10* 

w 

/(/tee 

5  x  10* 

p 

rad/  tec 

1  x  10s 

1 

rad/  tec 

1  x  105 

r 

rad/ tee 

1  x  10s 

i. 

inchet /tec 

1  x  10* 

ic 

inchet  /  tee 

1  x  10* 

i. 

inchet / tee 

1  x  10* 

U 

ft/iec* 

5  x  10* 

V 

/(/tee* 

5  x  10* 

w 

ft/tee* 

5  x  10* 

t. 

inchet 

1  x  10’ 

»c 

inchet 

1  x  10* 

_ h _ 

inchet 

1  x  10* 

Figure  4.7:  Reduced  Order  Hover  Compensator.  The  poorly  damped  modes 
(a  =  -.5  ±  j'2.03)  dominate  the  pitch  angle  response  as  Figure  4.9  shows. 
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outer  loop  gains  were  set  using  the  RSANDY  program.  This  approach  did  not 
work  well  as  the  convergence  was  very  slow.  The  alternative  was  to  set  the  gains 
intuitively.  This  approach  was  actually  quite  reasonable  if  one  recalls  the  control 
structure  of  Figure  4.3.  The  gains  were  set  by  determining  how  much  velocity 
would  be  reasonable  to  use  to  correct  a  given  position  error.  For  instance,  if  the 
aircraft  were  10  feet  from  the  desired  hover  point  and  a  pilot  would  be  willing  to 
use  5  ft/ sec  of  y  then  Kiff  would  be  .5  ///,*“;■  This  approach  worked  well  for 
the  y  and  z  axes  but  failed  for  the  x  axis.  For  this  axis,  a  conventional  transfer 
function  analysis  was  used  to  set  the  5  outer  loop  gains  which  were  then  adjusted 
in  simulation. 

4.5.1.  Transfer  Function  Analysis  for  X  Axis  Outer  Loop 

The  first  step  in  the  process  required  a  transfer  function  from  ucommand  to 
uacllI0i  for  the  helicopter  with  velocity  inner  controller.  This  came  from  use  of 
Bernard’s  code  for  calculating  the  matrix  of  transfer  functions  for  any  MIMO 
linear  dynamic  system.[7]  These  transfer  functions  were  11'*  order  so  to  make  the 
process  tractable,  the  NAVFIT  program  at  NASA  Ames  simplified  them  to  3rd 
order.  The  following  transfer  function  resulted: 

Ugttual  _  _ ;£! _  (4  4) 

Ucomm  and  (s  +  .24)  ^(s  +  .5)  +  1.872j 

To  justify  the  use  of  9  and  q  as  second  and  third  derivatives  of  x,  consider  the 
longitudinal  equation  of  motion  from  Etkin:  [14] 

Fj  -  mg  sin9  =  m  [u  +  (gf  +  q)  w  -  (rf  4-  r)  v]  (4.5) 

where  if  the  following  assumptions  are  made: 


RD-H151  946  A  MODERN  CONTROL  DESIGN  METHODOLOGY  WITH  APPLICATION  TO  2/3 
THE  CH-47  HELICOPTEROJ)  AIR  FORCE  INST  OF  TECH 
URIGHT-PATTERSON  AFB  OH  R  D  HOLDRIDGE  JAN  85 
UNCLASSIFIED  AFIT/CI/NR-85-33D  F/G  1/4  NL 
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<?§  *3  0  Earth  rotation  negligible 

rf  ss  0  Earth  rotation  negligible 

qw  fa  0  Second  order  effect 

tv  fts  0  Second  order  effect 

sm  0  «  9  Small  angle  assumption 

F„  sy  0  Reasonable  for  a  helicopter  in  hover 

ti  «  x  Small  angle  assumption(0) 

then  these  simplications  result: 

x  ss  —g9 

(4-6) 

u  «  —g9  S3  —gq 

The  actual  gain  setting  is  done  by  including  the  following  for  the  PIDD2D 3 
controller  in  Figure  4.8. 


P I DD2D 3  =  Kpx  +  Kf — (-  Kdxs  +  Kptxs2  +  Kp»xss 

3 


where 


Kp  =  Kitj 

K[  =  KiJj 
K0  =  Ka 
Kp  i  =  Ka 

Kd .  =  Kiq 


Rewriting  Equation  4.7  as: 


PIDD'D 1  =  Kn.  (V  +  +  £-)- 

\  An *  An*  An*  An*/  S 
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u-  forward  velocity 

Figure  4.8:  Transfer  Function  Analyaia.  The  PIDDiDi  compensation  was  cal¬ 
culated  to  cancel  the  two  lightly  damped  poles  of  Equation  4.4  and  move  the  two  poles  at 
the  origin  to  the  left. 

then  the  four  numerator  zeros  of  the  PIDD*D 3  controller  were  selected  to  cancel 
the  lightly  damped  poles  of  Equation  4.4  and  to  draw  the  two  poles  at  the  origin 
to  the  left.  The  gain  KDS  was  selected  for  good  speed  of  response.  Following  are 


the  gains  calculated: 


KP  =  KXp  =  -1.085 


K[  =  Kxi  =  -.131 


Kd  =  Kxd  =  -2.94^ 
Kd, i  =  Kx$  =  .551^ 


:g  fP a 

(ft  error ) 

fps 

(ft  sec  error) 
jp° 


(4.10) 


Kn*  =  Kxo  ~  .393 


(deg/ sec) 


4.5.2.  Outer  Loop  Simulation 


With  these  gains  as  a  starting  point,  the  time  responses  were  improved  using 
the  onboard  simulation  in  the  flight  computer.  Although  this  approach  (“tweek- 
ing”  the  gains)  may  seem  somewhat  unscientific,  it  was  appropriate  for  several 


reasons: 
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•  Extensive  use  of  the  onboard  simulation  had  the  added  advantage  of  helping 
discover  many  errors  in  the  flight  software  before  actual  flight. 

•  Pilot  comments  concerning  the  performance  and  response  characteristics 
could  be  better  incorporated  into  the  design. 

•  Working  directly  with  the  fixed  point  digital  computer  avoided  the  addi¬ 
tional  time  and  effort  required  to  digitize  and  scale  the  continuous  design. 

•  The  transient-free  switching  between  the  velocity  command  and  position 

hold  modes  could  be  developed.  This  is  discussed  in  the  next  section. 

The  PID  gains  coming  from  this  simulation  are  shown  below: 

*■*<■  =  --75  !7i^T  Krr  =  -*•«  !7l£=T  k»  =  -*•<>  Tfiu^S 

KX,  =  -3.8  x  10-  Ky,  =  -7.6  x  10-  |  Kz,  =  -1.9  x  10- 

Kxd  =  -3.0  to  Kyd  =  -.18  to  Kzd  _  -4.5  to 

«x»  =  .2  to  =  0.0  jg 

^<3  =  2.0  Krp  =  0.0 

(4.11) 

Time  responses  from  simulation  are  shown  in  Figures  4.9  to  4.13.  These  figures 
show  the  results  of  both  the  inner  and  outer  loop  designs. 

4.6.  Flight  Test  Implementation 

The  flight  implementation  of  the  hover  control  system  was  based  on  the  soft¬ 
ware  and  flight  procedures  developed  for  the  longitudinal  CAS  flight  test.  As 
before,  the  sensor  data  were  filtered  by  the  TR-48  analog  computer  before  being 
digitized  and  sent  to  the  Sperry  digital  computer.  There  were  two  areas  where 
the  hover  controller  was  quite  different  from  the  longitudinal  system  and  required 
new  flight  capabilities.  The  first  was  the  inertial  position  and  velocity  data  and 
the  second  was  the  transient-free  switching  required  to  make  the  transition  from 
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Figure  4.9:  Hover  Forward  Velocity  Step  Command 

poor  damping  is  evident  in  the  pitch  angle  response. 


in  Simulation.  The 


I  DEG)  I  INCHES)  (FEET/SEC) 

10.0  5.0  0.0  5.0  10.0  -l.o  -O.s  0.0  0.5  l  .0  -10. 0  -S.o  0.0  5.0  10.0 
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Figure  4.10:  Hover  Side  Velocity  Step  Command  in  Simulation.  S 

velocity  performance  is  adequate  but  the  poorly  damped  pitch  mode  is  also  excited. 


lULG)  l INCHES)  lEEET/SCa 

10. U  5.0  0.0  5.0  1U.0  1.0  -0.5  0.0  0.5  1.0  -10.0  -5.0  0.0  5.0  10.0 
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Figure  4.11:  Hover  Vertical  Velocity  Step  Command  in  Simulation.  The 

heave  response  is  adequate  and  does  not  excite  the  pitch  modes  as  strongly  as  the  lateral 
velocity  response  of  Figure  4.10 
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Figure  4.12:  Hover  Forward  Position  Step  Command  in  Simulation.  This 
response  came  from  the  onboard  simulation  and  shows  the  the  outer  loop  performance 
using  the  flight  software. 
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Figure  4.13:  Hover  Lateral  Position  Step  Command  in  Simulation.  Again, 
the  flight  software  shows  good  outer  loop  performance. 
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pilot  velocity  command  to  automatic  position  hold. 

4.6.1.  Inertial  Velocity  and  Position  Data 

Although  the  hover  controller  uses  the  inertial  data  (z,  z,  y,  y,  z,  and  z )  as  it 
does  the  other  measurements,  a  considerable  effort  was  required  to  get  these  data.1 
Since  the  INS  positions  drifted  so  quickly  and  there  were  no  inertial  velocities 
available  from  the  INS,  an  alternative  source  for  these  data  was  needed.  The 
ground  based  tracker  at  Crows  Landing  was  able  to  provide  these  data  using  a 
ground-to-air  telemetry  link  that  was  specially  developed  for  this  program.  The 
steps  required  to  make  these  ground  based  position  measurements  usable  by  the 
control  laws  are  described  below: 

•  The  laser  or  radar  tracker  measured  position  of  the  aircraft  in  a  runway 
based  polar  coordinate  system. 

•  These  measurements  of  azimuth  angle  [Az],  elevation  angle  {El),  and  range 
(r)  were  telemetered  from  the  ground  tracking  station  to  the  helicopter. 

•  These  data,  as  well  as  tracker  status  information,  were  decoded  and  scaled 
into  units  common  to  the  rest  of  the  flight  software. 

•  The  Az,  El,  and  range  data  were  converted  to  a  runway  based  rectangular 
coordinate  system  with  a  new  origin  located  over  the  runway. 

•  These  data  were  used  with  aircraft  accelerometer  data  (rotated  into  the 
runway  coordinate  frame)  in  a  second  order  complementary  filter  to  esti¬ 
mate  z,  x,  y,  y ,  z,  and  z.  Figure  4.14  shows  the  block  diagram  of  this 

_ complementary  filter. 

'This  work  was  done  primarily  by  Bill  Hindson  of  NASA  Ames 
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limiter 


Figure  4.14:  Second  Order  Complementary  Filter.  The  filter  uses  acceleration 
and  position  to  determine  the  smoothed  position  and  velocity  needed  by  the  control  law. 

•  These  smoothed  values  of  x,  x,  y,  y,  z,  and  i  were  then  rotated  through 
the  aircraft  heading  angle  to  the  heading-oriented  inertial  frame  required 
by  the  control  system. 

•  Based  on  the  tracker  status  information  coming  from  the  ground  and  based 
on  data  reasonableness  checks,  an  algorithm  kept  the  inertial  data  consis¬ 
tent  during  short  term  tracker  breaklocks.  For  longer  term  breaklocks,  the 
experimental  control  system  was  disengaged  to  avoid  the  large  control  mo¬ 
tions  caused  by  trying  to  follow  bad  data.  Initially,  the  laser  tacker  was 
used  since  it  provided  more  precise  range  information  (1-2  foot  acurracy). 
Unfortunately,  the  laser  had  frequent  and  unpredictable  breaklocks  which 
made  the  data  essentially  unusable  in  the  control  loop.  Because  of  this  in¬ 
ability  to  hold  lock,  the  radar  tracker  was  used  for  the  flight  test  although 
its  accuracy  was  only  5  -  10  feet. 
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4.6.2.  Transient-Free  Switching 

Since  this  control  system  had  both  manual  and  automatic  capability  in  the 
three  body  axes,  a  way  was  needed  to  transition  smoothly  among  these  different 
control  modes.  This  task  was  complicated  by  the  following  characteristics  of  the 
control  system: 

•  The  pilot  had  the  freedom  to  change  heading  at  any  time. 

•  The  x  coordinate  of  the  desired  hover  point  had  to  follow  the  i  coordinate 
of  the  actual  position  when  the  pilot  was  commanding  i  velocity.  At  the 
same  time,  the  helicopter  had  to  hold  both  y  and  z  position  in  the  heading 
inertial  frame. 

•  Same  as  above  in  the  y  and  z  directions. 

•  A  detent  on  the  pilot  controls  was  needed.  If  the  pilot’s  control  was  less 
than  the  detent  value,  that  axis  was  in  position  hold  mode;  else,  the  pilot 
was  commanding  a  velocity. 

Figures  4.15,  4.16,  and  4.17  show  the  switching  logic  for  the  three  axes.  The 
assembly  code,  shown  in  Appendix  I,  implements  this  logic. 

4.7.  Flight  Test  Results 

The  hover  flight  testing  was  done  at  the  Crows  Landing  test  facility.  The 
testing  was  limited  to  this  location  since  the  system  required  the  use  of  the  radar 
tracker  at  Crows.  Unlike  the  longitudinal  CAS  control  system,  the  hover  controller 
was  very  difficult  to  debug  and  make  operational.  The  flight  testing  was  divided 
into  three  phases  to  accommodate  these  difficulties.  The  first  phase  developed  the 
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Figure  4.15:  X  Axis  Transient- Free  Switching  Logic.  The  transformation  from 
body  to  runway  uses  heading  angle  from  the  INS.  The  best  deadzone  or  detent  value  was 
about  .25  inches. 
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Figure  4.25:  Hover  Forward  Velocity  Step  Command  in  Simulation.  The 

pitch  angle  response  is  improved  from  the  original  design  of  Figure  4.9. 
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quicken  the  lateral  response  while  in  hold  mode  hence  reducing  any  coupling  due 
to  action  in  the  longitudinal  axis.  With  these  two  changes  to  reduce  the  coupling, 
and  the  redesigned  inner  loop  velocity  system,  the  final  flight  testing  began. 


4.7.4.  Final  Closed  Loop  Flight  Test  in  Hover 

With  the  redesigned  controller,  the  velocity  performance  was  significantly  im¬ 
proved.  Figure  4.28  shows  the  response  to  a  command  in  x.  The  poor  pitch 
damping  has  been  eliminated  and  the  coupling  to  bank  angle  is  gone.  The  y 
axis  velocity  performance  remains  good  as  shown  in  Figure  4.29.  The  z  velocity 
response  remained  almost  identical  to  the  original  design  shown  in  Figure  4.21. 
The  position  hold  performance  is  also  evident  in  these  figures  when  the  velocity 
commands  are  removed  and  the  system  reenters  the  postion  hold  mode.  Figures 
4.30  and  4.31  confirm  the  good  hold  performance  in  the  y  and  z  axes  but  the  poor 
damping  in  position  hold  in  x.  The  x  axis  position  hold  dynamics  are  dominated 
by  a  slow,  poorly  damped  mode  (<;  w  .4  and  u  w  50  sec).  A  significant  amount 
of  flight  time  was  spent  adjusting  gains  in  the  PID  outer  loops  to  improve  this  x 
hold  performance.  Shown  below  are  the  final  set  of  outer  loop  gains  which  resulted 
from  these  efforts.  Later  flight  tests  used  the  integrator  in  the  x  axis  PIDD2D 3 
controller  only  when  the  error  was  less  than  40  feet.  This  improved  the  damping 
slightly  to  about  .5. 


Kxp  =  -  20 
KX,  =  -7.0  x  10-* 

Kxd = -2  0  ft 
Kx,  =  10  ft 
Kxq  -  3  0  (sfe 


Kyf  =  -1.0 


Kzp  —  —2.0 


( ft  error ) 


Kyi  =  -7.6  x  10  3  j/(  ,{?’„ror)  KZI  =  -1.9  x  10~3  ^ 
KYd  =  -10  ^  KZd  =  -4.5  ^ 

*»>  =  -■  5S 

Kyp  =  _-5  y 

(4.12) 


t  9tc  error ) 
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Figure  4.24:  Redesigned  Hover  Compensator.  Only  the  columns  of  K  associated 
q  and  9  and  the  row  of  C  corresponding  to  longitudinal  control  {Se)  are  changed  from  the 
initial  design  of  Figure  4.7. 
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good  time  responses.  In  this  case,  the  measurement  noise  characteristics  were  left 
unchanged  and  an  unrealistically  high  value  of  the  vertical  velocity  disturbance  was 
used.  Specifically,  the  vertical  gust  root  mean  square  (rms)  was  increased  from 
2.3  ft/ sec  to  10  ft/sec.  Vertical  gust  was  selected  since  it  affects  the  pitch  angle 
more  strongly  than  the  other  disturbances.  With  this  one  change,  the  RSANDY 
program  was  used  to  find  a  new  compensator.  To  speed  up  the  convergence  in 
the  RSANDY  program,  only  the  columns  of  the  Jfm,„  matrix  associated  with 
measurements  of  q  and  9,  and  the  row  of  Cm,n  associated  with  the  longitudinal 
control  were  allowed  to  vary.  This  approach  was  also  logical  since  we  wished  to 
keep  the  vertical  and  lateral  axes  unchanged  from  the  first  design.  The  redesigned 
compensator  in  shown  in  Figure  4.24  and  can  be  compared  to  the  initial  design  in 
Figure  4.7.  Figure  4.25  shows  the  simulation  response  of  the  redesigned  velocity 
command  inner  loop  with  the  improvement  in  pitch  damping  compared  to  the 
initial  design  shown  in  Figure  4.9.  Figures  4.26  and  4.27  show  that  the  y  and  z 
responses  were  essentially  unchanged  by  the  redesign. 

Since  there  was  nothing  in  the  simulation  to  suggest  that  there  would  be 
coupling  from  i  command  to  <P,  the  approach  to  solving  this  problem  was  based 
on  the  experience  gained  thus  far.  Two  changes  were  made  to  the  controller  which 
would  have  to  wait  for  flight  to  be  evaluated.  The  first  change  was  the  zeroing  of 
the  feedforward  gain  from  i  command  to  <5„  in  the  N  matrix  of  Figure  4.3.  This  was 
a  logical  approach  to  solving  the  problem  since  the  N  matrix  was  highly  dependent 
on  accurate  modeling  and  the  longitudinal  flight  test  had  already  shown  the  model 
to  be  lacking.  The  other  change  made  to  solve  this  coupling  was  to  include  nonzero 
values  for  Ky *  and  Ky  P  in  the  lateral  PID  outer  loop.  This  change  was  made  to 
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Figure  4.23:  Preliminary  Lateral  Position  Step  Command  in  Flight.  The 

y  position  hold  performance  is  well  damped  and  similar  to  the  simulation  of  Figure  4.13. 
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Figure  4.22:  Preliminary  Forward  Position  Step  Command  in  Flight.  The 

x  position  damping  is  very  poor  (f  w  .1),  unlike  the  near  critical  damping  of  the  simulation 
shown  in  Figure  4.12. 
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Figure  4.21:  Preliminary  z  Flight  Response.  The  5  second  time  to  steady  state 
matches  the  simulation  shown  in  Figure  4.11. 
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Figure  4.20.  Preliminary  y  Flight  Response.  The  y  response  is  well  behaved  and 
similar  to  the  simulation  results  of  Figure  4.10.  The  peak  roll  angle  is  about  8  degrees  for 
both  responses  to  a  command  of  about  10  ft) see. 
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Figure  4.19:  Preliminary  z  Flight  Response  with  Pitch  Roll  Coupling. 

With  the  instability  of  Figure  4.18  corrected,  the  roll  coupling,  shown  here,  was  discovered. 
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ary  i  Flight  Response  with  Pitch  Oscillation.  The 

q  response  of  Figure  4.9  became  unstable  in  flight  as  the  pitch 
ling  using  the  outer  loops  worked  well. 
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system  in  flight.  Figure  4.18  shows  the  flight  results  with  the  extremely  poor  x 
axis  performance.  After  a  redesign  of  the  velocity  inner  loop  system,  described  in 
the  next  section,  another  problem  was  found.  There  was  unacceptable  coupling 
between  the  x  and  y  axes.  The  coupling,  evident  in  Figure  4.19,  was  manifested 
as  roll  osscillations  resulting  from  the  x  command.  The  next  section  also  describes 
the  approach  taken  to  solve  this  pitch  to  roll  coupling  problem.  The  coupling 
was  one  way,  however,  as  seen  in  Figure  4.20  where  the  response  to  a  step  in  y 
is  quite  acceptable  and  similar  to  the  simulation  results  of  Figure  4.10.  The  z 
command  capability  is  also  quite  good  as  Figure  4.21  shows.  The  x  position  hold 
performance  of  Figure  4.22  was  very  poor  due  to  the  low  damping  («  .1).  Y 
position  (Figure  4.22)  was  much  better  damped  and  faster  than  x.  The  z  position 
hold  performance  was  very  good  with  vertical  position  changes  of  less  than  10  feet 
during  the  velocity  commands  of  Figures  4.19,  4.20,  and  4.21.  Use  of  the  radar 
tracker  data  in  the  inner  loops,  which  was  considered  risky  due  to  its  complexity, 
worked  well  throughout  the  flight  test. 

4.7.3.  Hover  Controller  Redesign 

The  redesign  of  the  system  was  necessitated  by  bad  performance  in  two  modes. 
First,  the  x  velocity  response  was  slightly  unstable  in  flight.  The  other  problem 
was  the  coupling  from  x  velocity  command  to  roll  angle.  The  first  problem  was 
handled  by  redesigning  the  inner  loop  velocity  control  system  in  order  to  slow  the 
longitudinal  response.  This  redesign  was  first  attempted  by  changing  the  weighting 
matrices  in  the  RSANDY  program  to  get  a  better  damped  longitudinal  response. 
This  approach  did  not  work  so  the  technique  of  Section  3.4.1  was  used.  Section 
3.4.1  described  using  an  arbitrary  set  of  measurement  spectral  densities  to  achieve 
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new  capabilities  needed  by  the  hover  controller  including  data  uplink  capability 
and  complementary  filtering  to  get  smooth  inertial  data.  The  second  phase  in¬ 
cluded  preliminary  flight  test  which  discovered  poor  velocity  performance  which 
necessitated  a  redesign  of  the  velocity  inner  loop.  The  final  phase  of  the  flying 
evaluated  the  redesigned  control  system. 

4.7.1.  Support  Systems  Development  Flying 

The  complexity  of  the  hover  controller  required  that  essentially  all  the  aircraft 
systems  and  all  the  ground  support  equipment  be  working  in  order  to  exercise  the 
system.  A  number  of  flights  was  required  just  to  ensure  that  the  uplink  system 
and  the  associated  complementary  filters  were  producing  good  inertial  data.  Once 
these  systems  were  operating  correctly,  the  flight  testing  continued  with  checks 
of  the  mode  switching  and  transient  suppression  logic  while  using  the  real  data 
coming  from  the  complementary  filters.  It  was  while  doing  this  work  that  the  laser 
tracker’s  poor  ability  to  hold  lock  was  discovered  and  the  decision  was  make  to  go 
with  the  less  accurate  radar  tracker. 

4.7.2.  Preliminary  Closed  Loop  Flight  Test  in  Hover 

Preliminary  closed  loop  testing  included  velocity  step  commands  in  the  three 
axes  and  changes  in  desired  position  while  remaining  in  the  hover  hold  mode. 
These  closed  loop  tests  confirmed  what  the  longitudinal  CAS  tests  had  already 
shown.  The  flight  responses  were  less  damped  than  the  simulations  had  predicted. 
In  other  words,  we  couldn’t  achieve  as  high  a  bandwidth  in  flight  as  in  simulation. 
This  was  most  evident  in  the  x  axis  where  the  well  damped  velocity  response 
in  the  simulation  (Figure  4.9)  turned  into  a  neutrally  stable  or  slightly  unstable 
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Figure  4.16:  Y  Axis  Transient-Free  Switching  Logic.  The  switching  logic  for 
the  x  and  y  axes  was  identical.  The  best  detent  was  .25  inches. 
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Figure  4.26:  Hover  Side  Velocity  Step  Command  in  Simulation.  The 

response  is  nearly  identical  to  the  original  design  of  Figure  4.10. 
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Figure  4.27:  Hover  Vertical  Velocity  Step  Command  in  Simulation.  The 

response  is  nearly  identical  to  the  original  design  of  Figure  4.11. 
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4.8.  Summary  of  Results  of  the  Hover  Controller 
Design 

As  with  the  longitudinal  CAS,  the  discussion  of  results  is  separated  into  two 
groups: 

•  the  effectiveness  of  the  methodology 

•  the  flight  test  results 

The  hover  controller  emphasized  the  usefulness  of  the  design  methodology  for  a 
more  complicated  control  system.  To  have  used  classical  incremental  loop  closures 
to  do  this  design  would  probably  have  taken  longer  or  would  have  required  more 
specific  experience  in  helicopter  control  systems  than  I  had.  This  task  also  showed 
the  advantage  of  using  a  modern  control  inner  loop  to  modify  the  open  loop  plant 
in  such  a  way  as  to  increase  the  physical  intuition  for  the  design  engineer.  The 
increased  physical  intuition  made  classically  designed  outer  loops  simpler.  In  this 
case,  the  plant  was  changed  from  control  motion  in,  measurement  out  to  desired 
output  in,  actual  output  out.  This  change  simplified  the  selection  of  the  outer  loop 
control  structure  and  made  the  outer  loop  gains  more  intuitive.  Figure  4.3  showed 
these  advantages.  This  task  also  emphasized  the  relative  speed  and  ease  with 
which  design  iterations  can  be  made  on  MIMO  systems.  When  the  first  design 
of  the  hover  controller  was  found  unacceptable  in  flight,  the  redesign  described 
in  Section  4.7.3  was  done  in  only  2-3  days  which  avoided  delays  in  the  flight 
testing.  As  with  the  longitudinal  CAS  design,  the  analysis  tools  developed  to  use 
the  methodology  (described  in  the  various  appendices)  were  sufficient  but  their 
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Figure  4.28:  Final  x  Flight  Response.  The  instability  in  pitch  (Figure  4. 18)  and 
the  pitch  roll  coupling  (Figure  4.19)  are  gone  but  the  pitch  angle  damping  is  still  less  than 
the  simulation  response  of  Figure  4.25. 
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Figure  4.29:  Final  y  Flight  Response.  The  y  response  is  little  changed  from  the 
first  flight  test  results  of  Figure  4.20. 
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Figure  4.30:  Final  Hover  Forward  Position  Step  Command  in  Flight.  The 

performance  is  much  improved  from  Figure  4.22  but  the  damping  of  ps  .4  is  still  not  good 
enough  for  operational  use. 
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Figure  4.31:  Hover  Lateral  Position  Step  Command  in  Flight.  The  lateral 

position  response  is  adequate  and  little  changed  from  the  first  tests  of  Figure  4.23. 
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“user  friendliness”  needed  improvement. 

In  flight  test,  the  hover  controller  was  fairly  successful.  Many  advanced  hel- 
copter  or  VTOL  airplane  designs  call  for  a  translational  velocity  control(TVC) 
system  such  as  was  tested  here.  Normally  the  evaluations  of  these  concepts  sel¬ 
dom  leave  the  simulators  to  address  the  hardware  and  software  difficulties  of  flight 
implementation.  This  test  reemphasized  two  of  the  important  difficulties  of  TVC 
systems,  the  inertial  position/velocity  sensor  and  the  human  factors  involved.  The 
primary  contribution  of  this  work  was  the  development  of  a  flexible  TVC  system 
where  these  type  issues  can  be  studied.  Specifically,  this  system  showed  good  ve¬ 
locity  command  performance  in  all  three  axes,  excellent  hold  performance  in  the 
lateral  and  vertical  axes,  and  marginally  acceptable  hold  performance  longitudi¬ 
nally.  The  switching  logic  worked  well  from  a  control  viewpoint  but  the  pilots  who 
0'  ’  flew  the  system  commented  on  the  need  for  a  better  indication  of  switching  from 

hold  mode  to  velocity  command  mode  in  each  axis. 


Chapter  5. 
Conclusions 


The  conclusions  are  separated  into  those  applicable  to  the  design  methodology 
and  those  associated  with  the  flight  b;sts. 

5.1.  Methodology 

•  The  process  of  scaling  and  using  the  modal  input/output  measures  is  an 
effective  way  to  reduce  the  order  of  the  compensator. 

•  The  scaled  block  minimal  realization  of  the  compensator  is  useful  in  iden¬ 
tifying  unimportant  measurements  and  controls. 

•  The  decoupling  feedforward  matrix  depends  heavily  on  an  accurate  model 
so  practical  designs  will  usually  require  some  sort  of  integral  control.  For 
very  poor  models,  integral  control  alone  should  be  used  for  implementing 
output  commands. 

•  The  software  tools  were  adequate  for  application  of  the  methodology. 

•  Both  “modern”  and  “classical”  control  techniques  are  important  for  MIMO 
control  system  design.  The  specific  application  determines  the  appropriate 
techniques  to  use.  In  this  research,  the  use  of  a  modern  control  inner  loop 
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wiih  classically  designed  outer  loops  was  a  useful  approach  for  the  hover 
controller. 


6.2.  Flight  Tests 

•  Although  requiring  some  change  of  pilot  technique  (retraining),  the  decou¬ 
pled  velocity  and  climb  rate  controller  was  well  received  by  the  pilots  who 
flew  it. 

•  The  hover  controller  performed  adequately  as  a  translational  velocity  com¬ 
mand  system,  had  good  position  hold  capability  in  the  vertical  and  lateral 
axes,  but  its  hold  performance  in  the  longitudinal  direction  was  marginal. 

•  Integral  control  was  crucial  to  achieve  decoupled  control  for  both  the  cruise 
and  hover  control  systems. 

5.3.  Lessons  Learned 

Finally,  two  “lessons  learned”  (or  relearned)  during  this  research  should  be 
emphasized,  even  though  they  may  seem  obvious.  First,  there  is  no  substitute  for 
experience.  For  this  methodology,  experience  was  important  in: 

•  selecting  the  correct  units  for  scaling  the  dynamic  system 

•  determining  which  gains  are  “small”  for  compensator  simplification 

•  selecting  outputs  and  their  weightings  in  the  optimization  using  the  ROPT- 
SYS  and  RSANDY  computer  programs 

•  selecting  scale  factors  for  fixed-point  scaling,  Appendix  E 
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•  selecting  stick  and  collective  lever  gains  for  the  pilot 

•  determining  the  structure  of  the  integral  control  loops 

Also,  experience  in  use  of  the  methodology  itself,  especially  the  design  tools,  was 
critical.  The  hover  controller,  though  much  more  complicated,  took  about  as  much 
time  to  design  as  the  longitudinal  CAS.  The  other  lesson  is  that  the  design  of  the 
control  logic  is  often  the  easiest  and  fastest  step  in  building  an  operational  control 
system.  Most  of  the  work  is  spent  on: 

•  software  design,  coding,  and  testing 

•  hardware  modifications  and  testing 

•  ground  based  closed  loop  testing 


0) 


Chapter  6. 

Recommendations  for  Further  Research 


The  CH-47  research  helicopter  at  Ames  is  a  very  flexible  test  vehicle  and  is 
being  improved  by  the  addition  of  a  floating  point  digital  computer  programmable 
in  a  higher  order  language.  With  this  improvement,  a  number  of  potential  research 
projects  should  be  considered: 

•  Parameter  identification  to  improve  the  existing  models  used  for  design  and 
simulation 

•  Refinement  of  the  two  designs  presented  here  and  pilot  evaluation  in  a  more 
realistic  setting  such  as  instrument  landing 

•  Outer  loop  guidance  work  (Microwave  Landing  System,  4-dimensional  nav¬ 
igation,  etc.)  using  these  inner  loops 

•  Application  of  singular  value  LQG-LTR  (Linear  Quadratic  Gaussian  -  Loop 
Transfer  Recovery)  to  account  for  unmodeled  rotor  dynamics 

One  difficulty  in  applying  the  methodology  was  the  poor  convergence  character¬ 
istics  of  the  first  order  gradient  algorithm  in  the  RSANDY  program.  A  second 
order  technique  to  speed  convergence  would  be  an  important  improvement  to  the 
program.  Another  possibility  for  research  is  finding  a  way  of  commanding  a  sys¬ 
tem  without  exciting  all  the  closed-loop  modes  similar  to  the  method  described  in 


120 


Cha.pter6 


121 


Chapter  3  for  the  case  of  the  full-order  compensator. 

Saberi  has  shown  a  technique  for  calculating  helicopter  stability  derivatives 
during  low  speed  flight  near  the  ground. [18]  This  research  vehicle  is  an  excellent 
testbed  for  validating  these  derivatives. 
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Appendix  D. 

Set  Point  Design 


This  appendix  derives  the  feedforward  matrix  which  commands  a  dynamic 
system,  including  compensator,  to  a  new  equilibrium.  This  matrix  turns  desired 
outputs  into  the  steady-state  controls  needed  to  achieve  the  outputs.  The  restric¬ 
tions  are  that  these  outputs  (or  new  operating  point)  be  physically  realizable  and 
that  the  number  of  controls  be  equal  to  the  number  of  outputs. 

Consider  the  following  dynamic  system: 

x  =  Fx  +  Gu 
y*  =  H,x  -f-  Z?,„u 

( D.l ) 

z  =  Az  +  By,  +  Gzue 
u  =  Cz  +  Dy,  +  ut 

where: 

x-  plant  states 
y,-  measurements 
z-  compensator  states 
u-  controls 
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0J' 


The  general  form  of  the  transformation  is: 

Ti  0  ...  O' 

0  T2  0  0 

:  0  0 

0  0  0-*- 

77*  o  ...  o- 

o  r2- 1  o  o 

:  0  0 
0  0  0  ci. 

where: 

cl.ii(fi9i  -  e*hi)  eiCt,2i  -  Mf,2t-1  1 

rp  _  .  C/.2«-i  (/«?«’  ^«^i)  ff»'^f.2i’-l  ^iCf.2i-l  J 

—  y»C/.2«  —  1 2  +  (e*  —  ^i)c/,2i-lC/,2»  +  <7,C[,2,2 


(C.ll) 


(C.12) 


t-  the  i"1  complex  mode 

/-  row  index  corresponding  to  the  largest  value  of  c\  +  c\  for  that  mode’s 
double  column,  or  the  largest  value  of  c\  for  the  a  real  mode 

A  listing  of  a  FORTRAN  subroutine,  MINCOM,  which  does  this  transforma¬ 


tion  is  shown  in  Appendix  L. 
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The  “0  1”  in  the  Cmin  matrix  results  from  scaling  the  system  by  the  largest  values 


in  the  C  double  columns. 

Consider  a  second  order  system  in  general  form: 

[tH;  flfcMI;]* 


u  =  ( cx  C2\  Zl 
.  z2  . 


(C.  6) 


The  desired  transformation  will  put  this  system  into  the  following  form: 


ai  a2J  l  z: 


u  —  f  0  1 


(C.  7) 


The  similarity  transformation  requires  the  two  systems  to  have  identical  eigenval¬ 


ues,  that  is: 


I  si  -  A\  =  |s/  —  i4min| 


ax  =  fg  —  eh 
a2  =  e  +  h 

Introducing  the  transformation  matrix  and  expanding: 

T Amin  -  AT 


(C.8) 


Cmtn  =  cr 


r„  r12|f  o  1  =  re  /]  rr„ 

T2 1  T22\[fg-eh  e  +  h\  L0  h\  Lt21 


I  0  1 1  =  c,  c2 


Tn  Ti2 
Tn  T22 


(C.9) 


where  cl  +  c]  has  the  largest  magnitude  of  any  row  pair  in  double  column  associated 
with  the  mode  being  made  minimal.  Solving  the  equations  above  we  have  the 


desired  transformation: 


c2{fg  -  eh)  ec2  -  6c, 
-Cjjfg-eh)  gc2-hci 
-fc]  +  (e  -  h)cic2  +  gc\ 


(C.  10) 


.  r  •"  -  «r  1 
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where: 


A  = 


★  ★ 
★  ★ 

0 

0 


0 

★  ★ 
★  ★ 


0 

0 


0 

0  0 

B  =  [Full] ,  r  x  p 

(7  =  [Fuff] ,  m  x  r 

We  want  a  transformation  to  a  new  form: 

Z  —  Amin  Z  +  B„ 


0 

★ 


,r  x  r 


n-C„ 


where: 

Amin-  minimal  form  of  A 
Bmin~  minimal  form  of  B 
Cmin  -  minimal  form  of  C 


Amin  — 


Bmin 


0  1 
a,  a* 

0 


*  * 
*  * 


0 

0  1 
Oj  02 


,n  x  p 


,n  x  n 


★  *  0  1 

0  1** 
*  *  *  * 


,  m  x  n 


z  —  Tz' 

Amin  =  T~lAT 
Bmin  =  r->£ 

cmin  =  cr 


(C.2) 


(C.3) 


» 


(C'A) 


(C.5) 
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Appendix  C. 

Minimal  Realizations 


The  design  methodology  described  in  Chapter  2  used  minimal  realizations 
in  two  places.  The  first  was  when  the  ROPTSYS  computer  program  displayed 
the  compensator  in  minimal  form  to  be  better  suited  for  the  optimization  in  the 
RSANDY  program.  This  eliminated  redundant  parameters  which  could  cause 
trouble  in  the  RSANDY  gradient  search  procedures.  The  second  use  of  a  mini¬ 
mal  realization  came  when  the  discrete  compensator  was  transformed  to  minimal 
form  for  computational  efficiency.  In  the  first  case,  the  transformation  was  from 
arbitrary  form  to  block  modal  form  then  to  block  minimal  form.  The  second  was 
from  an  arbitrary  2x2  block  form  to  the  block  minimal  form.  The  derivation  is 
shown  for  an  arbitrary  2x2  system  then  expanded  for  any  order. 


Given  the  following  form  of  the  dynamic  system: 


z  —  Az  +  By 


u  =  Cz 


(C.l) 
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By  the  duality  property  of  regulators  and  estimators  (Figure  B.l),  these  gains, 
KT  =  (R  4-  NQNt)~1(NQTt  +  HmP),  are  determined  using  the  randomly  dis¬ 
turbed  equations  of  motion: 


x  =  Ft  +  Tw 
ym  =  Hm  +  Nw  +  v 


(B- 12) 


where 


Q-  noise  spectral  density  matrix  of  plant  distrubances,  w 
R-  noise  spectral  density  matrix  of  measurement  noise,  v 


Regulator 

F  G  H  L  A  B  S  C 

Estimator 

FT  HmT  r T  NT  Q  R  P  KT 

Figure  B.l:  Duality  Between  Regulators  and  Estimators.  This  shows  the 
property  of  duality  which  allows  the  use  of  the  regulator  results  for  design  of  an  optimal 
estimator  (a  steady-state  Kalman  filter). 
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Adjoining  the  constraints  (the  equations  of  motion)  to  form  the  Hamiltonian: 


H  =  t  +  A  T(Fx  +  Gu) 


(B.  5) 


Recalling  the  optimality  conditions: 


A”  = 


(*«) 


Introducing  the  system  equations,  x  =  F x  +  Gu,  and  expanding  the  optimality 
equations,  we  have  the  Euler-Lagrange  equations  (here  in  matrix  form): 

xl  _  F  -  G{B  +  LT  AL)~l  LT  AM  —G(B  +  LT  AL)~lGT  )\x 

Aj  -  [—M7 AM  -  M7- AL(B  +  LTAL)'LTAM  -F  -  M7" AL(B  +  LTAL)~lGT  j  LA 

(H.8) 

As  shown  in  Bryson  and  Ho  [4]  or  Franklin  and  Powell  [6],  the  solution  to  these 
equations  is  A  =  Sx  where  S  =  A-Xll.  A_  and  X -  are  the  submatrices  of  the 
eigenvector  matrix  of  the  Hamiltonian  matrix  associated  with  eigenvalues  having 


negative  real  values,  i.e. 


[xi  rr.  r+m-i 

LaJ  -  La_  a+jU+J 


(B.9) 


With  this  solution  for  A,  the  optimal  steady  state  control,  u,  can  be  expressed  as 
a  linear  combination  of  the  state  variables,  x: 


u  =  Cx 


C  =  (B  +  LT  AL)~l{LT  AM  +  GT  S) 


(B.  10) 


The  same  approach  applies  to  finding  the  estimator  gains,  K,  of  the  equation: 


x  =  Fx  +  Gu  +  K[ym  -  Hmx) 


(B.  11) 
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Optimal  Compensator  Design 


The  design  methodology  described  in  Chapter  2  uses  an  optimal  full  order 
compensator  as  the  starting  point.  This  appendix  summarizes  the  derivation  of 
the  optimal  compensator.  From  Hall  and  Bryson,  we  see  that  to  design  a  set 
of  regulator  gains  which  minimize  a  quadratic  performance  index,  we  minimize 
the  Hamiltonian  with  respect  to  the  control. [17]  In  this  case,  the  performance 
index  includes  the  control  in  the  output,  thus  enabling  the  weighting  of  state 
rates(accelerations).  This  is  essential  in  aerospace  applications  where  vehicle  ac¬ 
celeration  is  an  important  parameter  in  the  design  and  analysis  of  the  control 
system.  Starting  with  the  modified  performance  index: 

J=  f  Cdt  =  f  \{yT  Ay  +  uT  Bu)d.t  (5.1) 

Jo  Jo  2 

where: 

y  =  Mi  4-  Lu  (5.2) 

yT  =  xTMT  +  uTLT  (5.3) 

Expanding  the  integrand  of  J: 

c=l-  [xtMtAMx  +  xTMTALn  4-  u T LT AMx  +  uT(LT AL  +  5)u]  (5.4) 
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< 7 -  standard  deviation  of  the  noise  variable 
TV-  correlation  time  of  the  noise 

With  these  transformations,  the  resulting  compensator  will  use  scaled  mea¬ 
surements  to  calculate  a  scaled  control  signal.  If  we  want  to  use  the  compensator 
in  the  physical  system,  we  need  only  unscale  the  gain  matrices.  The  compensator 
based  on  the  scaled  variables  is: 

i  =  Az  +  B$, 

-  *  (-A-®) 

G  =  Cz+  Dtf, 

We  unscale  the  system  by  replacing  scaled  vectors  G  and  jf,  with  their  unsealed 
equivalents  G  =  Te~lu  and  y,  —  Tm~ly 

i  =  Az+  B  Tm~l  y, 

(A- 9) 

G  =  TeCz  +  D  Tm~l  y. 

Now  the  compensator  uses  actual  (unsealed)  measurements  and  gives  unsealed 
control  signals  as  outputs. 

To  make  the  scaling  process  consistent,  I’ve  listed  some  rules  of  thumb  below: 

•  Scale  the  matrices  consistently;  for  example,  if  a  measurement  is  also  a 
state,  use  the  same  units. 

•  Scale  intermediate  state  variables  in  an  actuator  model  the  same  as  the 
control  itself.  For  example,  if  we  have  a  first  order  actuator  model,  ua  = 
-ati„  +  au„  then  scale  the  actuator  position  state(u„),  its  rate(u„),  and  the 
command(ut)  identically. 

•  Similarly,  if  sensor  noise  filters  are  included  in  the  plant  model,  then  these 
r  rise  filter  states  should  be  scaled  the  same  as  the  measurements  they  filter. 
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turbances: 

y,  =  Tmg, 

Vt  =  Tpye 

(AA) 

w  =  Tdw 
v  =  Tmx) 

With  these  transformations,  the  scaled  dynamic  system  is: 

x  —  Fi  +  Gu  +  f  u> 

y,  =  ff,£  +  D,uu  +  tftZ)  4-  0  (A.5) 

yc  =  He2  +  Deua 

where: 

F  =  T,~lFT, 
g  =  rrlGTe 
r  =  T.-lTTd 
a,  =  Tm-lH,T, 

D.*  =  Tm~lD,uTe 

(^.6) 

N  =  Tm~lNTd 

B(  = 

=  Tp-lDeuT' 

Q  =  Tf'QTf1 
R  =  Tm~lRTm~l 

The  scaled  power  spectral  density  matrices  were  derived  using  the  approxima¬ 
tion: 

PSD  »  2<r2Te  (A.7) 

where: 
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f  -  plant  dynamics  matrix,  n  x  n 
G-  control  distribution  matrix,  n  x  m 
r-  plant  disturbance  distribution  matrix,  n  x  m' 

H,-  state  to  measurement  distribution  matrix,  p  x  n 

Dtu-  control  to  measurement  distribution  matrix,  p  x  m 

N~  plant  disturbance  to  measurement  distribution  matrix,  p  x  m' 

He -  state  to  output  distribution  matrix,  p'  x  n 
Deu~  control  to  output  distribution  matrix,  p'  x  m 
Q -  plant  disturbance  spectral  density  matrix,  m'  x  m' 

R-  sensor  noise  spectral  density  matrix,  p  x  p 

The  scaling  process  continues  by  describing  the  changes  of  units  on  the  states, 
controls,  etc.  as  simple  transformations.  For  instance,  if  we  want  new  states,  i , 
and  new  controls,  u,  to  be  $*,12,  •••  3xnxn  and  s„,ui>  sUju2,  ...  sUnu„,  then 

we  can  define  scaling  (also  similarity)  transformations: 


x  =  T,x 
u  =  Teu 


{A.2) 


I 


A 

B 

A 


( 


(  £) 

L« 


Appendix  A. 

Engineering  Scaling 


This  appendix  derives  engineering  scaling  equations  used  in  the  ROPTSYS 
computer  program.  This  process  transforms  the  model  of  the  physical  system  into 
a  “similar”  model  where  the  units  of  the  variables  have  changed.  Similar  means 
the  eigenvalues  of  the  system  are  not  changed  by  the  transformation  to  the  new 
coordinates.  As  described  in  section  2.2,  the  new  units  are  chosen  to  make  the 
new  variables  of  the  dynamic  system  of  equal  importance  to  the  design  engineer. 
The  process  begins  with  the  linear  model  shown  below: 

x  =  Fx  4-  Gu  4-  Tw 
y,  —  H.x  +  D,uu  4-  Nw  4-  v 

(A.l) 

ye  =  Htx  +  Deuu 

J  =  f  {yj  Aye  +  uTBu)  dt 

Jo 

where: 

x-  system  states,  nx  1 
z-  compensator  states,  r  x  1 
u-  controls,  m  x  1 
w-  plant  disturbances,  m1  x  1 
y,~  sensor  measurements,  p  x  1 
ye~  weighted  outputs,  p'x  1 
v-  sensor  noise,  p  x  1 
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If  DD, 0  =  0,  then  these  equations  can  be  rewritten  as: 


1  fzl 

G 

JU  + 

.  BDIU  +  Gt . 

F  +  GDH,  GC 

3(1  +  DtuD)H,  A  +  BD, 

Defining  xt  =  ,  we  can  rewrite  the  equations  above  as 


U' 


(D.  2) 


Xt  —  FjXt  4-  G  rue 


(D.3) 


Expressing  the  desired  outputs,  t/£>,  as  a  linear  combination  of  i r  and  u,  we  have: 


j/d  =  Hdxt  +  Ld^c 


(DA) 


At  steady  state,  x V  =  0,  and  the  two  previous  equations  become: 


Ft 

Gt 

XT  _ 

0 

Hd 

Fo . 

.ujf, 

Vd 

inverting: 


it 

Ft 

Gt 

-i  j 

’  0  ' 

.  U'  . 

1  Ho 

Lp . 

■  Vd. 

*  Af 

0  ' 

*  N . 

.Vd. 

(D.5) 


(D.6) 


The  steady-state  controls  are  uf  =  Nyo  and  new  equilibrium  state  vector  is  xT  = 


Myo  where: 


(0.7) 

This  last  equation  is  used  by  the  SETPNT  program,  Appendix  L,  to  calculate  the 
N  matrix. 


Appendix  E. 

Fixed  Point  Scaling 


This  appendix  describes  the  technique  of  scaling  the  analytical  designs  to  run 
on  the  Sperry  18 19 A  flight  computer.  The  process  is  similar  in  principle  to  scaling 
engineering  problems  for  an  analog  computer.  This  computer  is  an  18  bit  fixed- 
point  digital  computer.  There  are  two  problems  which  must  be  considered  when 
doing  this  scaling.  The  first  is  avoiding  overflows  (exceeding  217  —  1  during  cal¬ 
culations)  and  the  second  is  maintaining  precision  in  the  results.  The  procedure 
which  follows  handles  both  these  potential  problems. 


Consider  the  compensator  dynamic  system: 

z  =  Az  +  By 
u  =  Cz 


(F.l) 


where  y  and  u  are  in  engineering  units  (not  yet  computer  scaled).  In  the  computer, 
these  variable  have  computer  scaling  factors,  K,  such  that  zK ,,  zKt,  yKv,  and 
uKu  have  units  of  bits.  For  example,  if  Kt  =  500  then  then  5  deg  of  9  is 
2500  bits  in  the  computer.  In  these  computer  scaled  variables,  the  compensator 
appears: 


!*!  K,  =  Ka  [X|  -  M  K,  +  Kb  |B|  [y|  K, 


k,ka 

M  K„  =  K„  [O]  W  K, 


(E.2) 
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where  A,  B,  and  C  are  in  engineering  units.  The  scale  factors,  K}  for  y  and  u 
are  part  of  the  computer  environment  (set  up  by  the  programmers  of  the  original 
flight  program)  but  we  need  to  calculate  K,  and  Kt.  This  is  simplified  by  the 
digital  implemention;  we  have  z*+ 1  and  z*  rather  than  z  and  z.  This  means  only 
one  scale  facter,  Kt,  is  needed.  To  find  this  factor,  first  estimate  the  largest  value 
that  any  z*+i  or  z*  can  achieve  by  finding  the  maximum  single  product  in  the 
matrix  multiply,  [5]  [y]m0J.  For  controls,  measurements,  and  desired  outputs,  the 
maximum  values  can  be  set  using  engineering  judgement  and  intuition.  Since  we 
also  want  precision  in  the  z  term,  we  select  Kt  so  that  zmaM  uses  all  of  the  18  bits 
available: 

217 

K,  = -  (E.  3) 

Zmai 

where  zmaj  is  rounded  up  to  the  next  power  of  2.  217  is  used  since  the  largest 
negative  number  expressed  in  18  bits  is  -217. 

With  the  Kt  term,  the  problem  of  overflow  is  solved.  Now  we  need  only  ensure 
that  precision  is  maintained  in  the  calculations  by  choosing  the  additional  scale 
factors,  KA,  KB,  and  Kc  that  scale  the  elements  of  the  A,  B,  and  C  matrices. 
Making  use  of  all  18  bits,  we  can  find  these  scale  factors  in  the  same  way  as  the 
Kt  term  above  was  calculated: 


>17 


tf.A  = 

kb  = 

Kc  = 


1  mn*  K, 

917 


L  K, 


(EA) 


>17 


~'matKu 


where  amaj,  bmaj,  and  cm0,  are  the  maximum  elements  of  the  A,  B,  and  C  matrices 


which  have  been  rounded  up  to  the  nearest  power  of  2. 
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One  product  from  the  matrix  multiplies  is: 

(E.s) 

This  number  is  included  in  the  double  precesion(36  bits)  A  register  in  the  18 19 A 
and  is  always  less  than  235  -  1.  We  want  to  accumulate  these  double  precision 
elements  to  get  one  element  of  By.  Finally,  we  divide  by  KA  (or  equivalently  shift 
the  A  register)  to  regain  the  single  precision  inner  product.  This  is  done  for  each 
of  the  matrix  multiplies. 

The  same  approach  is  used  to  scale  the  feedforward  matrix  and  the  integral 
gain  matrix  used  in  the  longitudinal  CAS.  The  two  BASIC  computer  programs, 
which  do  this  scaling  for  the  longitudinal  CAS  and  for  the  hover  controller,  are 
listed  in  Appendix  N  along  with  example  data  files. 


Appendix  F. 


CH-47  Research  Helicopter 

The  helicopter  used  for  this  research  is  a  highly  modified  version  of  the  Boeing* 
Vertol  CH-47  “Chinook”  used  by  the  U.S  Army  for  cargo  and  troop  transport. 
Figure  F.l  shows  the  tandem  rotor  helicopter,  which  is  operated  by  the  NASA 
Ames  Research  Center.  Reference  10  is  a  more  complete  description  of  this  partic¬ 
ular  helicopter  including  the  many  modifications  made  to  the  basic  CH-47.  Below 

are  listed  some  of  the  modifications  and  improved  capabilities: 

•  Full  authority,  variable  stability,  fly-by-wire  flight  control  system  in  all  four  axes. 

•  Programmable  analog  and  digital  computers  capable  of  executing  the  control  laws. 

•  Programmable  force-feel  system  on  the  experimental  pilot’s  stick. 

•  Flight  instrumentation  system  capable  of  recording  over  100  variables  at  100  times 
per  second. 

•  Operator’s  console  for  control  of  the  experimental  systems. 

•  Additional  sensors:  INS,  radar  altimeter,  body-mounted  accelerometers,  improved 
air  data  sensors,  numerous  control  position  sensors,  boom-mounted  angle  of  attack 
and  sideslip  vanes,  rate  gyros 

•  Digital  ground  to  air  uplink  capability 

Figures  F.2  and  F.3,  from  reference  10,  show  the  cabin  layout  in  this  experi¬ 
mental  vehicle  and  a  block  diagram  of  the  experimental  control  system. 
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Figure  F.2:  Cabin  Layout.  The  research  helicopter  requires  a  crew  of  4;  safety  pilot, 
experimental  pilot,  research  system  operator,  and  crew  chief. 


Figure  F.3:  Experimental  Flight  Systems.  The  flexibility  of  the  experimental 
system  is  emphasized  here  where  we  note  the  many  interfaces  between  the  various  com¬ 
ponents. 
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Appendix  G. 


CH-47  Linear  Models 


The  models  described  in  this  appendix  were  calculated  using  the  information 
from  reference  11.  The  general  forms  for  the  decoupled  4th  order  models  and  for 
the  coupled  order  are  shown  in  Figures  G.l  and  G.2.  The  linear  models  for 
the  flight  conditions  related  to  this  research  are  shown  in  Figures  G.3  to  G.16. 


Figure  G.l:  Longitudinal  and  Lateral  4th  Order  Models.  The  longitudinal 
model  above  was  used  in  the  design  of  the  longitudinal  CAS  (Chapter  3). 
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Figure  G.3:  Linear  Model  for  Airspeed  of  00  knots,  Climb  Rate  0  ft/min. 
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Figure  G.5:  Linear  Model  for  Airspeed  of  80  knots,  Climb  Rate  0  ft/min. 
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Figure  G.6:  Linear  Model  for  Airspeed  of  60  knots^  Climb  Rate  500 
ft /min. 


150 


Figure  H.2:  C 

importance  of  t 


Figure  H.l:  Bessel  Filter  Analog  Flow  Diagrams.  These  filters 
grammed  on  the  airborne  TR-48  analog  computer. 
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shown  in  Figure  H.2  where  the  unfiltered  and  filtered  measurements  are  compared. 
By  the  end  of  flight  test,  all  these  filters  were  replaced  by  hardwired  3rd  order  Bessel 
filters  located  in  a  signal  conditioning  box. 


Appendix  H. 
Bessel  Filters 


The  Bessel  filters  described  here  were  designed  to  eliminate  the  “3  per  rev” 
and  “6  per  rev”  harmonics  at  11  Hz  and  22  Hz  due  to  the  225  rpm  rotor.  The 
break  frequency  was  chosen  at  5  Hz  as  a  compromise  between  noise  attenuation 
and  measurement  bandwidth.  The  actual  filter  designs  came  from  reference  15. 
Nine  filters  were  patched  on  the  airborne  TR-48  analog  computer: 

4th  order-  Body  axis  accelerations  (.4*,  A¥,  At) 

3rd  order-  Roll  rate  (p) 

2nd  order-  Pitch  rate  ( q ),  Yaw  rate  (r),  Velocity  (u),  Altitude  ( h ) 

The  transfer  functions  for  these  filters  are  shown  below: 

974603 

s<  +  98ss  4-  4323s2  4-  96906s  4-  974603 
_ 29220 _ 

s3  4-  76s2  +  2374s  4-  29220  (tf.l) 

987 

s2  4-  54s  4-  987 

Figure  H.l  shows  the  analog  patch  diagrams  for  these  filters.  Their  effectiveness  is 
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Figure  G.15:  Linear  Model  for  Airspeed  of  20  knots,  Climb  Rate  500 
ft/min. 
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Figure  G.14:  Linear  Model  for  Airspeed  of  20  knots,  Climb  Rate  0  ft /min. 
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Figure  G.8:  Linear  Model  for  Airspeed  of -20  knots,  Climb  Rate  0  ft/min. 
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Figure  G.7:  Linear  Model  for  Airspeed  of  00  knots,  Climb  Rate  -500 
ft/min. 
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Appendix  I. 
Flight  Software 


This  appendix  lists  the  flight  computer  software  unique  to  this  project.  It  is 
written  in  Sperry  1819A  assembly  code  and  internally  documented.  Three  different 
sets  of  the  code  are  shown.  The  initialization  code  was  executed  whenever  the 
experimental  control  system  was  disengaged.  It  did  such  things  as: 

•  reset  the  compensator  states  to  zero 

•  set  the  trim  values  of  controls  and  states  to  the  current  values  of  these 
variables  in  preparation  for  system  engage 

The  experimental  controller  subroutines  computed  the  control  laws,  described  in 
Chapters  3  and  4,  and  were  executed  at  20  Hz.  Both  the  longitudinal  CAS  and 
hover  controller  are  shown.  The  last  item  shown  in  this  appendix  are  the  in¬ 
structions  which  reserve  memory  for  the  variables  and  contants  unique  to  the 
subroutines. 
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Initialisation  Code,  (contd) 
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•  i 


-*  b  0.  J  b  «Q  1 

Ink  o  1 

b  n  N  b  n  Uft 

V-  J  b  O  fcl 

a.  *  «  *  o.  K  ob 


1  n  aft  ■*  k  a  h 

: *  l *  5f  sif t  P  ^ 


.  gH  l  .*  l 

;  P-  |  P-!  ft 

i  «  x  U  »  w  x  u 

H  K  I  O  H  )C 

‘_G2  r  p  BS  * 


I  O  |  <0  «J 


H  *  o  <  a  ***  '  o  f  <«  P  o  «  a  ,  o> 

*  •  n  n  Mb  **  x  k  ■  n  n  a  ■  n  b  o  !  ■  b  t*  ft  o  [  ■  b 

tai  •-  oc  o  «•  w  1  n  n  o  n  -*  k  d  •*  y  h  h«P  •*  W  H  oU 

xb~b«*-.xoN  i  echsrmExh  b  ■  n  a  ••  o  k  ii  -*  k  — 

n  U  3  O-nn  W  x  ,  ua  I  ••  n  h  k  •«  r  ♦  is  •  o  P  m  p  ♦ 

H  *  pi  *■*  tft  P  n  *i  <-%  '  M  H  IQ  n  a  —  n  k«r*  —  ft*  fc  «  P  b  — 

w  u  -*  b  «-•  »-<  m  o  »m  in  n  u  n  b  n  u  «•  <k 

•  b  N  •  fc»  «j  —  K  l  *|«  K  H||  ulM  n  k  ^  b  -  b  o 

u  u  3  o  u  a*  a.  M  3  cxouUuuou  ♦pt*-*a  sob^anb  l  o  u  1 1  ^ 

kaziN^Lk  t  vimkkk  ih**E|  b  u  m  H  —  u  m  b  j  ■  n  «■*  a.  ■ 

I  kl«  «  fS  »  fc*M  M  ~  *■  h  ll  •«  M  M 

LNh^iiKClUN  NM*.Nr*«-«liatRbbO  -  (t  ».  B  S  I-  -<  k  k  ft  ft 

Z  K  N  N  L(  a  y  ft  3  KAHN  w  (ft  WftNN  UK  M  *  N  N  U 

lft.t-UOnfcto.Ni3  u3UZfc.NOUbZy.ND  *0  U  D  o  to.  N  O  ft  O  *  b  to.  n 

—  «  w  N  w  ^  ^  ^NwN««U*^wwM  u  b  **^WM  u  3  M  w.  w  H 

*  :  f  •  % '  % 


6-"-^ 
to  (O  fc  ft 
p  N  X  fcl 
K  «  t 


fcl  o  o  u 

toabo 

«K  U  U  U 
MOOD 
H  X  X  Z 
u  x  c  x 


-  -«g 

n  n  n  «■*  r*  ft* 

ft-  —  fc»  b  N  fe 

U  b  It  K  B  *  H 


O  *0 

O  X 

o 


-L  pn 


'(WIN 

O  ft* 

—  *  O  N 


(i»  b  N  it  fc  N  j  tfc  ft  ft*  tal  b  to.  •*  z  I  I  in  a. 

at  b  «  *  H  M  p  »  k  ft  i  k N  l  t K  b  b  N  b  at  b 

NNft  n  |  =>  (ft  H  p  N  ft  H  N  *  ft  ft*  N  N  b 

p  *  *  *  *m**  oxb  tcoK  i  i  kNNoaoo  okay 

e  I  :  ?  i 


N  o  f*  *  o 

y  *  •*♦<  a 

N  to*  N  to  U 

opofiJ 
ft  ft  a  u  u 
N  fc.  N  fc*  * 


jx  b  I  1  mb  I  x  1 1 

as  a  ct  u  u  (ft  a  ac  a  a  b  V  aft  D  a  i  u  ;  at  ft  b 

a  n  n  mPn  ,  a  H  n  n  n  a  j  n  « a x n 

at  t  ■  ft  at  ft  h  h  ft  ft  t  ft  ft  k  ape*  ebb  *  w  -  ►*  b  a.  at 

w  :  a.  *  a-  »-  a  a  b  b  b  n  o  ampmo  b  *  b  b  -  ►*  a  b  x  i  m  n  a 

a  n  aw  a  «o  a  o  H  a  a  a  y  tjafga  l  yj  a  a  9  a  a 

ea  '  o  i 

XI  i  |  O 

re  I  ••  if*  N  i  '  *  i  j 

n  bit*  P  -  a 

•to  5  N  N  o  i  "•  !  I  o 

b  b  b  jb  ID'  ft  !  « 

N  I  N  IN  < 

«  »  *  j  I  >  I  I 

•  L  «.  . ;  •!  •  I  !  -i  . 

^  JP»»*  *•  ••  o  >  mopar>  b»  mo*ab  4* 

2  y»  i  o  o  **  -  aapm  §  a  aoi*b»b*  p 

^  ^  J;  n  !  now  *  f  okn  *  fc  io  r»b  r-  < 

;  ;s-2:  ;  sisi  csssi  5  ;s:ss  ?  i 

«  f«**"  f-  f 

a  PPPf  p  b»  tn  b  m  »  •«#«  i  5  fc 

n  t  ttfct  #*•  p  n  •>,  r*  i  ■»  m  |kti»fcfc  •» 

r*  « 5?  ”  5  ?5W  '  prftpfftfc  h  <m  tnp  p 

o  o  o  o  o  o  o  o  c  o  >  o  o  o  o  >  o  >  e  >  o  a  a 

o  o  b  o  O  O  e  0*00  SoOoo  *  oooo*  f  ( 

m  r  wt  m  m  *  o  >Niaim4r<n  ft  m  *  o  o  n  .  «  m  *  trs 1  9  m  »  o  *  ••  m  m  « 

0  0  0*0*0  >000*0  (PftoSft  a«£«ltftl|ft2rif  fll 

»'r*#*»*  »•*  i  »•**  *ii  *•**•*  ^  #* 


uUuuu 

->  M  *J  J 

S6SS5 

*  a  \  k  n 

mrr 


H 


ft  «  N  W  N 

b  b  b  b  a 
ooo*o 
y  f  «f>« 


! 

♦  r- 


Appendixl 


169 


Longitudinal  CAS  Subroutine. 
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Longitudinal  CAS  Subroutine,  (contd) 
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Longitudinal  CAS  Subroutine,  (contd) 
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Longitudinal  CAS  Subroutine,  (contd) 
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Example  Data  File  for  Using  the  ROPTSYS  Computer  Program. 


:  ATA  FOR  OPTSYS  FOR  THE  CHAPTER  2  NAVION  EXAMPLE 

:  l;  inq.ir  iomf  iss.  im  itfi.dm.  ite3.  ifdfw.  ie.  idstb.  idbc.  iset.  ilnc  iscle  imin 

!  -'.0.0.  0.0.0  O  3. 0.0.0  2.0. 1.1.1 
NS  NC  NM  NPD  NO.NMOD.IREC 
■1  2  3  2.2.00 
:s  MATRIX (NS) 

1  0  1  O  01  01 

TC  MATP I  X (NC) 

01  10 

TO  MATR  I  X  (NPD) 

I  0  1  o 
TM  MATR  I  X  (NM) 

1010.  01 
TP  MATR  IX  (NO) 

10  10 
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ROPTSYS  User’s  Manual,  (contd) 


FM  rHT  REGULATOR  SYWTHFS I S  USE 
r  HD  FM  HS  AR  C  B 
F'-W  FILTER  SYNTHESIS  USE 
f  HO  CAR  Q  R 

DSTAB  R  FWTFR  DIAGONAL  ELEMENTS  ONLY 

(AJT  Vt|  fl  Q  ARE  DIAGONAL  OR  HILL  ( SELECT  TR) 

MO  IS  A  VECTOR 

IF  'AST  CARD  HX  l  IS  •  THTN  CONTINUE  TO  NEXT  CASE 

»*'-N  '  ^irp*cr>  flt!«  1  1  n#»  nptlan*  him  «rwl  *l'rl 


’'HROUTIME  LISTING 

SETUP  CAN  SET  UR  OPEN  LOOT  DYNAMICS  KART  I J 

INNER  OPTIMAL  CONTROLLER /OPT INAL  FSTIKATCR  DESTCN 
OPTIMAL  MODEL  FOLLOWING  DESIGN  OPTION 
CDtV-CCMPLEI  DIVISION 

AA/’RINT  matrix  printing 

again  reformats  eigenvalues  calculates  woal  surma  trices 

MULT  MATRIX  MULTIPLICATION 

madd  matrix  addition 
KIRN  matrix  transposition 

mjNV  MATRIX  INVERSION 

S'  ov  SOLVES  LYAPUNOV  EQUATION  FOR  SS  COVARIANCE 
***>€  ' OMPUTES  MCDAL  COORDINATE  MATRICES 
HER  FOR.HATS  C IGENVALUE/E  ICENVECTTR  OUTPUT 
TT  'TP^./TYS  TRANSFER  FlAICTtONS 
CHECK  CHECK  CONSISTENCY  Of  INPUTS 
POLES  PRINT  CUT  TT  POl.rS 
ZFRQ5  COMPUTE  SISO  TT  ZEROS 

ACOMR  COUVTT  ZEROS  BY  HROCXTTT '  S  METHOD  (  i  SUWCUTINTS) 

CCCRR 

SCL  PJNCTION 

res  id  f'vmm  itime  response)  Rrsrcurs 

AAJLANC  FIND  EIGENVALUES  AND  EIGENVECTORS  VIA 
>jp THE S  ■  THE  v*»  ALGCRITTM  (6  SUBROUTINES) 

■RERAN 

HQR2 

aalmak 

ScOH  CALCULATES  THE  MODAL  TORN  Of  THE  COWEMSATCR  («JUR  ItX^> 

MI  NOT*.  CALCULATES  THE  MINIMAL  FORM  Of  THE  OM>FNSATOR  (HDtLR  !  DCE) 
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Appendix  J. 

ROPTSYS  Computer  Program 


This  appendix  includes  the  users  manual  for  the  ROPTSYS  computer  program, 
which  is  located  on  the  FSD  VAX  at  NASA  Ames  Research  center,  and  an  example 
data  file.  The  data  file  was  used  to  calculate  the  full  order  compensator  in  the 
Navion  example  in  Chapter  2.  The  lengthy  FORTRAN  listing  is  not  shown. 
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Memory  Allocation,  (contd) 
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Hover  Controller  Snbrontine.  (contd) 
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Appendix  K. 

RSANDY  Computer  Program 


The  RSANDY  computer  program  is  a  modified  version  of  the  SANDY  com¬ 
puter  program  written  by  Uy-Loi  Ly  as  part  of  his  PhD  dissertation  at  Stanford 
[3]  and  later  modified  first  by  Gardner  [16].  It  is  stored  on  the  FSD  VAX  at  NASA 
Ames  Research  Center.  Ly  also  wrote  a  user’s  guide  for  the  SANDY  program. [8] 
This  appendix  gives  the  input  format  changes  to  make  the  SANDY  user’s  guide 
correct  for  the  RSANDY  program.  The  major  capabilities  added  by  the  RSANDY 
program  are: 

•  an  optional  gradient  step-size  reducer  from  Gardner  [16] 

•  a  linear  discrete  model  of  the  closed  loop  system  can  be  created  for  later 
simulation  studies 

•  a  leading  free  line  and  free  lines  before  all  the  data  items  are  included  to 
help  in  documenting  the  data  files 

The  changes  which  follow  apply  to  page  numbers  in  the  SANDY  User’s  Guide. 
Other  than  these  changes,  the  program  is  exactly  like  the  SANDY  program.  Also 
included  here  is  an  example  data  set  for  running  the  program.  This  data  was 
used  to  find  the  reduced  order  Navion  compensator  in  Chapter  2.  In  the  following 
changes,  the  new  variables  needed  by  the  program  are  italicized. 
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Change  1,  page  117a.  Running  the  RSANDY  Program 
@  RSANDY  Infile  Out  file  Simfile 
where: 

RSANDY-  A  VMS  command  file  which  runs  the  RSANDY.EXE  file  with  Infile, 
Outfile,  and  Simfile  as  data  files. 

Infile-  A  file  containing  the  input  data;  an  example  is  shown  at  the  end  of 
this  appendix. 

Outfile-  A  filename  where  the  program  will  write  the  output. 

Simfile-  A  filename  written  by  program,  if  IPLOT=l,  which  contains  the  lin¬ 
ear  simulation  models  used  by  the  SIMPLOT  program,  described  in 
Appendix  M. 

Change  2,  page  119a.  Item  1 
N  p ,  n ,  m ,  m ’ ,  p ,  p ’ ,  p ” ,  r ,  fl  ag ,  NNS,  IPPSS,  ICF,  ISS 
where: 

NNS-  Set  to  0. 

IPPSS-  Set  to  0. 

ICF-  Set  to  0. 

ISS-  Set  to  0. 

Change  3,  page  119a.  Item  2 

Maxfn.Nvar, To!, MSTEPMmearM Print, IDPRNJCLPRNMAPRNJPLOT,  DTJBG 
where: 

MSTEP-  Maximum  step  size  for  the  gradient  algorithm.  Start  at  100  and  make 
it  smaller  if  there  are  convergence  problems. 
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1DPRN-  =  0  for  no  input  data  printout. 

ICLPRN-  =  0  for  no  closed  loop  data  printout. 

MAPRN-  =  0  for  no  modal  analysis  printout. 

IP  LOT-  =  0  for  not  creating  a  simulation  model.  =  N  for  creating  a  simulation 
model  of  the  Ntr>  plant  condition. 

DT-  Cycle  time  for  the  discrete  simulation  model.  Rule  of  Thumb,  DT  « 

.2(_U_). 

Change  4,  page  126.  Add  Item  9 
Data:  XO 
Description: 

XO-  A  vector  with  (n  +  r  +  m')  zeros. 
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RSANDY  Example  Data  Pile. 


DATA  FOR  RSANDY  FOR  THE  CHAP  2  EX  THE  NAVI  ON  AT  tOO  KNOTS  (REDUCED  ORDER  CC**>) 
tPLNT  N.M.M*  P.P*  P"  R  FLAG  ftS  IPPSS  ICE  ISS 

Jvir  .Tol  KSTEpSNLrNEAR.  Tf°princ  IOPRN  ICLPRN  MAPRN.  IPLDT.  DT  IBC 
00  8  1  06  01.100.0.999  DO.  -1.  O  O.  O.  O.  Ob.  I 

PROBABILITY  DENSITY  OF  PLANT  CCWITICM  1 


1  O 

f  MATS IX  (N  X  N)U  W  Q.  THETA 

CHS  036  0.  32  2 
•  37 .2  02  176  0.0 
00191 .  -  0396.  -  2  96  0 
0.0. 1  0.0 
c  MATRIX (  n  x  m  ) 

O  1  0 
20  2.0 
•11  0.0 


0.0 

CAf*W  MATRIX (  n  X  •'  ) 

095  •  036.0.0 
37  2  02.0.0 
-  00191.  0396.0.0 
0.0  0.0 

HS  MATRIX  (  p  X  n  )  VELOCITY  HX7T 
l  0.0. 0.0 
0.1  0  0  .  176  O 
DSU  MATRIX  (  p  X  a  ) 

0.0 

o  o 

DSN  MATRIX (  p  X  •'  ) 

O  O  1  0.0 
O  0  O.  1  O 

HC  MATRIX  (  p'  *  n  )  VELOCITY  HX)T 

l  O  0.0.0 

O  - 1  O  0 .  1 76  O 

DCU  MATRIX (  p'  *  *  ) 

O  0 
O  0 

DTW  MATRIX (  p'  «  a  ) 

0  0  0  0  0 
0  0  0  0.0 

HP  MATR IX(  p '  '  k  n  1U.NQ  THETA  HDOT 

1  O  0  0  0 
O  1  0.0.0 
0  0  57  10 

O  O  O  57  3 

0  1  O  0.176  0 

DPU  MATR I X (  p' *  *  9  ) 

0  O 
0  0 
0  0 
0  0 
O  O 

DPR  MATRIX (  p-  X  a  > 

O  0  o.o 
O  O  0  o 
0  0  0  0 
0  O  0  0 

o  o  o.o 

uw  NOISE  TYPE  CRDER  RMS 

2  12  3 

FILTER  CONSTANT 


WW^MOISE  TYPE  ORDER  RMS 

2  l  2  1 


FILTER  CONSTANT 

U  MEAS  NOISE  TYPE  ORDER  RMS 
2  110 

FILTER  CONSTANT ( IKZ) 

FCOT  MEAS  NOISE  TYPE  CRDCR  RMS 

2  1.10 

FILTER  CONSTANT 
6  28 
ALPHA 
O  DO 

g  MATRIX (  p*  X  p'  ) VELOCITY. HX)T 

1  O  O 

0  1  O 

R  MATRIX (  a 
lOOOO  0.0 
O.  1  0 
A  MATRIX 
0  00006*00 
B  MATRIX 
O  1366 

2  676 
C  MATRIX 

0  00006*00 
13  14 

IDVAA (  Nver  ) 

3  4  5.6.7  0.11.12 

ITEM  9  XO 

0  0.0. 0.0. 0.0. 0.0.0  0  O  0.0 


i  X  a  )  ELEVATOR  .  THK7TTLE 


1  000 

0  2124 
3  241 

-0  1000E  01 
1  220 


13  59 


Appendix  L. 

SETPNT  Computer  Program 


This  appendix  lists  the  computer  program  which  calculates  the  feedforward 
matrix  described  in  Appendix  D.  This  version  does  not  have  the  Gt  matrix.  The 
data  formats  and  data  sequence  are  described  at  the  beginning  of  the  program. 
Also  included  is  an  example  data  set  (the  full-order  Navion  of  Chapter  2). 
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SETPNT  Computer  Program. 
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SETPNT  Computer  Program,  (contd) 
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SETPNT  Computer  Program,  (contd) 
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SETPNT  Computer  Program,  (contd) 
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SETPNT  Computer  Program,  (contd) 
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|.  I )  q,  I  !••»«  ((  ((  1 )b|pq 


call  aub 1  ( a 
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SETPNT  Computer  Program,  (contd) 
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SETPNT  Computer  Program  Example  Data  Set. 


N  M.P  R,  SAMPLE  TIME  EOR  SPHOV5  DAT 
4  2.3.4.  5 

F  MATR I X (  n  X  n  )- -U. W. Q , THETA 

-  04S, .036.0. -32  2 

-  37. -2  02. 176  0.0 
00191. - .0396. -2  98.0 

0,0, 1.0.0 

G  MATRIX (  n  X  m  ) 

O.  1 

-28  2.0 
-11  0,0 

0.0 

HS  MATRIX (  p  X  n  )- -U. HDOT, THETA 
1  0.0, 0,0 
O. -1  0.0. 176  0 
0.0.0. 1  o 

DSU  MATRIX.  PXM 

0.0 

0.0 

0.0 

A  MATRIX 

O  OOOOE  +00  1  000  0 .  000OE*0O  O  OOOOE-OO 

-14  69  -6  072  0  OOOOE *00  0  OOOOE ‘OO 

O  OOOOE *00  O  OOOOE  *00  -11  34  O  OOOOE  *00 

O  OOOOE *00  O  OOOOE'OO  0  OOOOE *00  -1  816 

B  MATRIX 

-02781  -O  8863E-01  6515 

1  485  0  8725  -33  66 

■0  1334E-01  -  2  6 22  2  441 

0.7711  -0  4641E-01  -10  10 

C  MATRIX 

0  OOO  -O  1000E-01  -0.  lOOOE'Ol  0  360OE-O2 
-O  1978  -O  2163E-01  0.3446E-01  -1  000 

D  MATRIX.  MXP 

0,0.0 

0  0,0 

HD  MATRIX.  M  X  N*R  -U .  HDOT.  THETA 

1 .0.0. 0.0. 0.0. 0,0 

0.  -  1  0.0. 176  0.0. 0.0.0 

LD  MATRIX.  M  X  M 

0.0 

0.0 
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SCALEM2  Computer  Program,  (contd) 
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SCALEM1  Computer  Program,  (contd) 


Z  Z  Z 

Ml- 

*  at  a.  £ 


£ 

£ 


-  i 


5  '  i  *  * 


* 

22- 


.  i 


,  i 

i/l  to. 


_  i 

t  z 

5! 


a  i  £ 
;  *  = 
su-u- 

K  •  K  Z 


£  e  v,t 

;£.£-£ 

v»  — 


cat 

aaE 


*  ^  g£  iu 


5ifee  e£  pp  -g  ??£ 
£  < 


r-  r  » 


cS 

3£ 

•  i 

|£ 

t 


— “  e  *  d  *  e*  *  es* 

-  pp£pi  pi  pi  pi  pfi 

~'Z  '  '  s'  S  s  S  '  S  * 

-*235“  cTSc*  :i  o*  £2 


„  _  z  -  _ 

11^5^  5 

2  S  *  z  <  <  <  * 


ss£ 


SeESE'S"  S’  E"  S'  EEE 

I. i  i£££i!I*_£i_£5_£5-£££ 

5Ji::£Sl*5  2Ji^i>8>--2i-S$t9S£S 

-»  —  --  *.  £  X  « 


e» 
pi 
*  ^ 

c* 

le 

£" 


el  el 
as  as 

ti  u 

le  ie 

f; 

o* 


a>  -,« 


*3£  - 


-  *- 


►p-p^ 

iifl- 

*ajjk 

mil 


-£  -£  jE 

=s  =s  A 
*'  P<  is 
a  2k  zi 

it  it 

=-  i- 


b« 

Pi 


cl  c  *  cc’io 


—  t--  H- 


s*  |g 


jjpjE 


a 


i  £ 


gfc 

1/3  Z 

UfA 

3*t/*r*H 

<  «  u  2 


5  £ 

£ 

/ 

« 

*> 

S’ 

—) 

s 

u.  >»  ® 

Spi 

M  S 

pi 

•»  s 
P3 

P~ 

!? 

P 

i 

a 

®  ® 

3  -3  ^ 

25  i 

s 

5 

H 

5 

t- 

5 

1“ 

3 

i*S 

23 

23 

?3 

f 

2 

2 

5 

| 

*-  -  t-  •* 

e  Mill 


aw 

«  ®  ®  - 


HV-at 


5553 

<<<< 

HHI/IW 


*7 


e»  o*  e« 

Pi  Pi  Pi 

c£  c*  o* 

|e 


S&5S& 


H  R  H  K  -> 

ppppI 

iiii{ 

<suZa 


j£  j£  ggggg 


s' 


E-  E-  ESS-E- 


-eS-r?: 

1*11*1) 


sss's! 
f*?S| SgB 

s,<sS;su 


g§§g*sg*s^Si?§SS2 


-  a 

h ;, 

s!  c: 
:£  25 

?:  ii 

—**•«*« 

S*$?£3 


34 

<5 

5  a5 


—C 
1-  £ 
;a 
£? 


3! 

!i 

?  z 

2$ 

it 


i 

ts 


y 

j=-S 

*S— T*.  — » 

.— t  fl  tb-  - - 

Si*  53 
3s:  |= 

2*c  S3 

Slc  ^ 

--*g 


i 

!i 


3  P 

c 

-8 


a 

K 

B 

| 

8 

E 

J? 


8? 

=e 


»e  p 

?L  if  if 


? 
s 

~e  e 
*S  i 


i  i' 


I! 


:  % 

1  ? 

2  P 


P 


pE  ? 

n  3  1/1 

—  i/>  £ 

<  -•£ 
5  S  w  k 


Si  '-i 

— ::ri  8§§|"5  =3  "3 


*»£S 


_  ^3  u"« z w £ S  13  a 

eE^;eeors;ee^s£-,;s;j5ur*8  si  ? 

._is?.-bp5._zs5:;s53sStt8?  rf  .  F 

■"  :-a  *  cisilli---  -- 

xiiizSissEuss 


poo 

lili 


?”  s*5 


s*2 


8  , 
s-  s-ig5  2 

?s  *?8*:r  3  £  : 

'y'jJ1  "*  ~T 

i  « 

"*^1 ''fei  ^  1 

M  iuC  PC  M  ~  “ 

i  ?e  s3  l  *5 

.1  Jl.lc  |  pi  «^s 

EE  EEEEe3„,“  lit  t“« 

eeueeEe^sSS*  s5?  s“* 

c»Sii28828SS9iS*?82SSg9 

•v'  .•  /  72«  444-C4-C-C  4  —  -  -  -*  “ 


(I  J|  <1S|I  i)  •  (2  W-W.ft2) 


SLALfMI  . . . .  t»)i>  t  tMNT  TAB|IO) 


AppendixN 


229 


SCALEM1  Computer  Program. 


n  & 


i  y  i  it 

•  s5  «»u  • 


=• 

-  aa$  —  lj 

*!ssU  I 

bbb^  3 

O  O  W  Jlu 

i- 


p-t;-e*^*g*P*ifi"s*e*  ;  BoBS-^S  s 
-B-B-B-S-g-S-S-B-S  >  u. 

ii!illi?isiiiiiiig=s|lis  s  s 

?8SSS!8!S188SSSt28SifIsSlllilll 


c  i 

u ;  —  w> 

-  ua  X 

8  o  !  §  6e  !  s  6  i  . 

I  m  §  £  *“  i  f  i  @2  s 

I  is-  *  g  g®  s  §  s  i  I 


2s  S  £  *“!  I  S  Sd5  P  -  P  t%S 

=5  i=  £  2=5  j  5  a:  a  -  a  *-Ss 

pt-  C<  m  x r  c  w  u _  —  ^ 

|*|  1“  *  ?  2*5  I  e  Si  SafesSsisS^ 

*£§  SB  siljJcSS  °i* 

ufS  ;£  £B«SS«=  eSs=«SSv,|£Sv,|5s3« 

£J  g|s  is  JUeSJh  as|fIs|i|i|i|BgS» 

£•  Ja  g|  B§25i8S-!8£S£Si£S8v.!S| 

s  sj.S  £1  e!!=:j*  ;s°35"5?is325“s4 

3s  pi§  =8  el’^ShsfiglighMg^ 


2  2-<  J  a  -£  ~  -  -  ¥-Sp 

(/>  X I  3  w  u  U ^  W 


c  c  ca  m 

- .  sut  e| 

-3  3t|  5* 
a”  “££  |g 
;g  BSq  s* 

2,  e==  -a 

as  b2.  <  . 

C5  ^;s  23 

:l  f* 

",  *¥S  ** 


iiBfcPP; 


»B  i°-  *5  e=  i  i  a  8j-  --fr  5 

“  =  *S£  »g!  u  3  3  3  hi  |  a  £  I  a 

sa  Sr:  32:  *  :  :  :  sfa  r  a  s  z  2 

N  !JS  y  ‘  1111*3  I 

sa  ggj*  gSs  *  .  .  u  «;2  £  s  :  s  S 


::::i:;:;:,“*‘*"***“**“*““****** 


ESSES! 

|  S  c  c  s  j 


9  Sfg 

uJ  at  i  i 

; 

£«  « ls£  * 

liisgii 


imi 

£11 

K‘P 

ri«M 


-  -gy£a>- 
8  8U~; 

:s«ib?5 

i:;# 

HIeeeee 

3 “-ococ a 
liSSS ?J?J 


Appendix  N. 

SCALEM1  and  SCALEM2  Computer 
Programs 


This  appendix  lists  the  computer  programs  for  doing  the  fixed  point  scaling 
described  in  Appendix  E.  Also  included  are  two  example  sets  of  data.  These 
programs  are  written  in  VAX  BASIC.  The  SCALEMl  computer  program  does  the 
scaling  for  the  longitudinal  CAS  controller.  The  SCALEM2  computer  program 
does  the  scaling  for  the  hover  controller.  The  description  of  required  data  is  at  the 
beginning  of  each  program.  The  SCALEMl  data  is  for  the  6'*  order  compensator 
with  integral  control.  The  SCALEM2  data  is  for  the  final  hover  controller. 
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S1MPLOT  Example  Data  File. 


4,8,0, .1. 15,0.0,0 
18,22,20.21.16, 17, 18.22 

'NAVION  WITH  FULL  ORDER  COMP  HDOT  COM  =  10  FT/SEC 
'FORWARD  VEL (FT/SEC)  ' 

'CLIMB  RATE (FT/SEC) 

'  PITCH  RATE (DEC/SEC) ' 

'  THETA  (DEG) 

'  ELEVATOR (DECREES)  ' 

'  THROTTLE (FT/SEC** 2)  ’ 

'FORWARD  VEL (FT/SEC)  * 

'CLIMB  RATE (FT/SEC)  ' 

0,0. 0.0,0  0, 10  0.0, 0,0. 0,0, 0,0. 0,0, 0.0. 0.0. 0.0. 0.0 
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SIMPLOT  Computer  Program,  (contd) 


DO  3100  1M.1WTS  _ 

IF  (yARAY  ( I)  LT  YWIN)  YMIN=YARAY  ( I ) 

V  SJl'pHTSOR  (JSORCINCRVE)  .YXWCIWCRVE)) 

CALL  AREA 2D (6  0  2  I) 

CALL  WA21E  I  'TIME  (.Pc)  ■ 

CALL  YVAWEitREriPLTIKMOm)  20) 

SS  SS 10  ,  SCALE'  TIME (HPtCTS)  W*  SCALE’  THAI) 
CALL  CURVE (TIME  Y AAA YNPNtS  O) 

CALL  ENDCR  (MPCHT) 

>400  CONTINUE 

CALL  FNCPL  (NPCNT) 

>r-00  CONTIWJE 

GO  TO  5000 

IMS  IS  THE  PART  THAT  DOES  4  PLOTS/FACE 

vnoo  DO  4500  NPOfT-1  MPACt 
CALL  PACE (11  0.8  5) 

^^NPLTR^i'oi'cALL  t**OT (  K3VIE  ' ) 

CALL  NO8R0R 

CALL  AREA 20 19  0.6  O) 

CALL  HEADlMfWEF  (CRT1T)  SO  2.1) 

CALL  E>CCR(NPCKn 

DJOU  START  THE  LOOP  FOR  PLOTTING  THE  4  CT^VES 
OO  4400  NCRVt=t  .4 

■Km  IS  WE  COUR*  or  THE  VECTOR  IN  DAT  SOW  plotted 

NCHT*  (4»NPCNT  4)  *MCRVE 
YHAX  -OAT  f  MCNT  *1.1) 

YWIN-YMAi 
DO  4100  l*l.MPNTS 

irc5ttl«tT*S?TaU  ’thar.yaratj.) 

IF  (YARAY(l)  LT  YM1N)  YM1N«YARM  (l» 

"ir0  S!I'?HT*»(««ORC((CRVE).Y«*CjNCRVE)) 

CALL  AREA 20(4  1.2  1). 

ssi  sasteffiV«)am).2o, 

CALL  GRAMO  .  SCALE'  TINE  (NEWTS)  TWIN  SCALE'  YMAX) 
C^  wSW  (TIME  .  YARAY  HwrfS.O) 

CALL  EHDCR(NPCWT) 

44CO  COKTltWE  _ _ 

CALL  EttJPL  (MPCHT) 

4500  CONTI** 

00  TO  5000 
MTO  CALL  DCA4EPL 
RETURN 
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SIMPLOT  Computer  Program,  (contd) 
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SIMPLOT  Computer  Program,  (contd) 
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m-  Number  of  controls. 

p"  -  Number  of  performance  variables. 

Item  3 

GRTIT,PLTIT(NPLT) 

where: 

GRTIT-  A  50  character  variable  containing  the  title  which  will  be  put  on  all  NPLT/NTYPE 
pages  of  plots. 

PLTIT-  An  array  of  NPLT  20  character  variables  containing  the  y-axis  labels  for  the 
plots.  The  x-axis  is  always  time. 

Note:  These  character  variables  must  be  exactly  50  and  20  characters  long. 

Item  4 

XC0(n  +  r  -f  ns) 

XCO  is  an  array  containing  n  +  r  +  ns  initial  conditions  on  the  simulation  states  in  the 
following  order: 


State 

Parameter 

H 

l 

plant  states 

In 

Zl 

l 

compensator  states 

ZT 

Wl 

l 

noise  filter  states 

AppendixM 


Description  of  “Plotdata”  File  inputs 


Item  1 

NT  YPE,  NPLT,  NPLTR,  D  Tl,  TMAX, IDF,  TR 

where: 

NTYPE-  “1”  for  1  plot  per  page,  horizontal.  “3”  for  3  plots  per  page,  vertical.  “4”  for  4 
plots  per  page,  horizontal. 

NPLT -  The  total  number  of  plots  to  be  made  NPLT/NTYPE  must  be  an  integer  less 
than  13. 

NPLTR  -  “1”  for  immediate  results  on  a  VT125  screen.  “0”  to  create  plot  files,  VECTR1.PLV 
and  PARM.PLV. 

DTI-  The  data  point  interval,  DTI/ DT(hom  RSANDY)  must  be  an  integer. 

TMAX-  Final  time  of  the  simulation,  TMAX/  DTI  must  be  less  than  5000. 

IDF-  “1”  to  rim  the  simulation  TR  seconds  prior  to  recording  data  for  plotting.  Thi3 
is  used  to  allow  filtered  noise  to  reach  a  steady  covariance.  “0”  to  start  taking 
plot  data  at  time=  0. 

TR-  The  number  of  seconds  to  wait  prior  to  recording  and  plotting  data. 

Item  2 

Plot  parameter  ID’s,  the  array  NPARM 

NPARM'is  an  array  containing  NPLT  identification  numbers  of  the  variables  to  be  plotted. 

The  parameters  are  identified  in  the  following  sequence: 


ID  number 

Parameter 

1 

1 

|  plant  states 

n 

Xn 

n  +  1 

Xl 

1 

1  compensator  states 

n  +  r 

Zr 

n  +  r  +  1 

Wi 

l 

1  noise  filter  states 

n  +  r  +  ns 

®r.i 

n  +  r  +  ns  +  l 

{ 

1  controls 

n  +  r  +  ns  +  m 

n  +  r  +  ns  +  m  +  l 

Vi 

i 

|  performance  variables(from  RSANDY) 

n  +  r  +  ns  +  m  +  p" 

Vp" 

where: 

n-  Order  of  the  plant. 
t-  Order  of  the  compensator. 

ns-  Number  of  states  in  all  the 

;  noise  filters. 
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Running  SIMPLOT 

©SIMPLOT  Plotdata  Simfile 
where: 

SIMPLOT-  A  VMS  command  file  which  runs  the  SIMPLOT.EXE  file  with  Plot- 
data  and  Simfile  as  data  files. 

Plotdata-  A  users  created  file  containing  data  used  by  SIMPLOT.  The  input 
item  descriptions  are  shown  on  following  pages. 

Simfile-  The  simulation  model  file  previously  created  by  an  RSANDY  run. 

The  SIMPLOT  program  will  ask  3  questions: 

1.  “Are  the  disturbances  random?” 

Type  “1”  for  yes  or  “0”  for  no.  If  no  noise  wa3  modeled  in  the  RSANDY  run  which 
created  the  Simfile ,  you  must  type  “0”.  If  you  had  noise  in  the  RSANDY  run  but  want 
a  clean  time  response,  you  can  type  “Cf  and  no  noise  will  show  up  in  the  time  responses. 
If  you  answered  “0”  to  this  question,  the  simulation  will  be  run.  If  you  answered  “1”,  the 
program  will  ask: 

2.  “Input  new  random  number  generator  seed?” 

Type  in  any  odd  integer  for  a  new  seed  or  type  in  zero  for  the  default  seed. 

3.  “Input  scale  factors  for  disturbances” 

To  force  the  simulation  to  be  driven  by  noise  having  the  same  RMS  as  the  RSANDY  data, 
type  “1,1,1.  •  1”  where  there  are  as  many  l's  as  there  were  random  inputs  in  the  RSANDY 
data,  i.e.  m' .  If  you  want  to  scale  the  RMS  of  the  noise,  change  the  “1  s”  accordingly.  For 
example,  typing  “2, .5"  will  make  the  first  noise  source  twice  the  RMS  of  the  RSANDY 
data  and  will  make  the  second  source  half  as  large. 


Appendix  M. 

SIMPLOT  Computer  Program 


This  appendix  is  the  user’s  guide  for  the  SIMPLOT  computer  program.  The 
simulation  part  of  the  code  was  written  by  Bruce  Gardner.  It  uses  models  created 
by  the  RSANDY  program  to  run  simulations  of  the  closed  loop  designs.  The 
program  has  the  option  of  including  the  random  disturbances  used  in  the  RSANDY 
design.  Up  to  12  variables  can  be  plotted  versus  time  and  displayed  on  a  VT125 
screen  or  sent  to  a  file  for  later  plotting.  The  listing  of  the  program  and  an  example 
data  set  (Navion  with  full-order  compensator  from  Chapter  2)  are  shown  at  the 
end  of  the  appendix.  The  user’s  guide  begins  on  the  next  page. 
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SCALEM2  Computer  Program,  (contd) 


r  *T  I 

2  24 O  PPtNT 

250  FOR  I- 
MT  I 

::6o  print 

22*0  PRINT 
2  280  PRINT 
'290  FOP  I 
NEXT  I 
>oo  PH  I  MT 
2  310  FCP  I 
NEXT  I 

2J2o  print 

2330  FCP  t 
serf  i 
2340  PatNT 
2350  PRINT 

2  360  pa  int 

2  370  FOP  1  = 
2  380  PRINT 
2390  FOP  l 

MT\  NEXT  I 

2400  print 

.'410  FOP  1  = 
NT\  NEXT  I 
2420  PRINT 
2430  FOP  1= 
NT\  NEXT  1 
2440  PRINT 
2450  FCP  I» 
2460  PP I MT 
2470  FCP  Ii 

247J  PPIVT 
2474  FCP  ti 
2480  PP I MT 
24<»0  PRINT 
2500  PPIKT 
2510  PPIKT 
2520  PPIKT 
2530  PPIKT 
2540  PPIKT 
2550  PPIKT 
2560  PPIKT 
2570  PRINT 
2580  PRINT 
2582  PPIKT 
2584  PRINT 
2586  PRINT 
2590  END 


TAB ( 10)  .  "THE  H  NATP 1 1  IS" 

1  TO  NC\  FCP  J=l  TO  MC\  PPIKT  TAB(J»8)  F>P(N(I  J)  >  \  NEXT  A  PPIKTv  K 


TAB ( 101  "THE  aXVUTT*  SCALED  VARIABLE  KATICES  APE" 

TAB  ( lO)  “THE  COHM/TFR  SCALED  8  MATRIX  IS" 

1  TO  FCP  J-l  TO  PPIKT  TABfJ‘8)  FUMBCSfl  J)  I  \  KFTT  A  PRINTN 
TAB  (lO).  "THE  CCAPVTCT  SCALED  C  MATRIX  IS" 

l  TO  MC\  FCP  J*1  TO  W\  PRIKT  TAB(J«0)  F>« (CCS  (t  J) )  \  NEXT  A  PRINT\ 
TAB (10)  "THE  OOfVTt*  SCALED  M  MATRIX  IS" 

1  TO  MC\  FCP  J*l  TO  NC\  PPIKT  TAB  (J*8)  .  F3«  (MCS  ( I  J) )  \  NEXT  A  PRINT\ 


TAB ( 10)  . 1 

TAB  ( lO)  .  ‘ 
1  TO  1*\ 

TAB  (10) 

1  TO  MC 


‘THE  COWLETELT  SCALED  MATRICES  APE" 

'THE  AF INAL  MATRIX  IS" 

PPIKT  TAB(  I  *0).  FTP  (AFIMAL  ( I ) )  .  \  NEXT  I 
‘THE  BF INAL  MATRIX  IS" 

FCP  J*l  TO  NHV  PRINT  TAB(J*8)  n«  (BFIMAL  (t .  J))  \  KEXT  A  PR  I 
’THE  CF INAL  MATRIX  IS" 

rrm  J»i  TO  f«\  PRIKT  TAB  (J*8)F>P(CF  INAL  (IJ))\  NEXT  A  PR  I 
THE  KFINAL  MATRIX  IS" 

FOP  J=1  TO  MC\  PPIKT  TAB(J*8)  .  FTP  (KFINAL  (I  J) )  \  NEXT  A  PR  I 


THE  KXFINAL  MATRIX  IS" 

PPIKT  TAB (t  *8)  .  FNP  (KXFTNAL  ( I )  )  ,  \  NEXT  I 
THE  OIFINAL  MATRIX  IS" 

PRINT  TAB(I  ♦§)  .  FFP  IKXIF INAL( I ) )  .  \  NETT  I 
THE  KXDf INAL  MATRIX  IS" 

PRINT  TAB  (I  *8)  .FRP  (KXDETMAL(l)  )  NEXT  T 


MC\ 

TAB (10)  .  ' 

1  TO  NC\ 

TAB  MO)  .  ' 

-1  TO  MC\ 

TAB  ( 10)  .* 

■1  TO  NC\ 

TAB  (10)  .  • 

■\  TO  NCy 

THE  A  MATRIX  PEQU 

THE  8  MATRIX  REQUIRES  A  LPTA' " ; 18 • SCALE 82  "THEM  STRAU' 

THE  C  MATRIX  REQUIRES  A  LPTA' 10  SCALEC2 . "THEN  STRAU" 

THE  N  MATRIX  REQUIRES  A  LPTA’  "  .  18  SCALE* J  ;  "THEM  STRAIT 

“THE  Oil)  ELEMENT  REQUIRES  A  LPTA' "  .  18-SCALEO  - 

THE  0(2)  ELEMENT  REQUIRES  A  LPTA '  "  .  18-SCALEO 

(3)  ELEMENT  REQUIRES  A  LPTA  ’ " .  18-  SCALEKX  (3)  .  "THEN  STRAU" 

‘  1  ELEMENT  REQUIRES  A  LPTA'"  10  SCALEOI  1)  .  THEN  STRAIT 

2  ELE«KT  REQUIRES  A  LPTA’"  18  SCALEOI  2)  "THEN  STRAU" 

3  CLEMENT  REQUIRES  A  LPTA'  1 8- SCALE O I  3)  "THEN  STRAU" 

1  ELEMENT  REQUIRES  A  LPTA IB-SCALEKXD  l)  .  “THEN  STRAU" 

2  ELEMENT  REQUIRES  A  LPTA' "  IB-SCALEKXD  2)  .  "THEN  STRAU" 

ELEMENT  REQUIRES  A  LRTA'"  10- SCALE KXD  3)  "THEN  STRAU" 


"THE  O 
"THE  KXi 
"THE  O! 
"THE  OI 
"THE  KXD 
"THE  ICO 
"THE  KXD 


31(1)  "THEM  STRAIT 
3(2)  :  "THEN  STRAU" 
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SCALEM1  Computer  Program  Example  Data. 


i) 


"NR.NM.NC.TP" 

6.4,2,  .05 

"A  MATRIX  (NR)" 

- . 23031.  89831. -  23642. 1  17826,  43492.  88649 
"B  MATRIX (NR  X  NM)  MEASUREMENTS  DISTRIBUTION  MATRIX  " 

-16  358. -456 .87. - . 2653.  05456 
-8. 274, -230.9, -. 15366.  012246 
10.224.195.99, .02073. -  005186 
2  39S48, 49  65, .00599.  000616 
2  749. 160  16. . 17429.  00936 
1  04027. - .01363. -  00002,  00428 

"C  MATRIX  (NC  X  NR)  COMPENSATOR  STATE  TO  CONTROL  DISTRIBUTION  MATRIX 
00.10.00.1.0,10.  688186 

-  91396. 2  70164. -  90986. 4  64211.  56013.1  0 

“N  MATRIX  (NC  X  NC)  COMMAND  DISTRIBUTION  MATRIX" 

-  08369. - .037896 

-  164,  07336 

"KI  MATRIX  (NC  X  NC)  INTEGRAL  CAIN  MATRIX" 

-  0042511. - .049046 

-  0040024. - .042269 


"BSCALE(NM) - THE  SCALINC  FACTORS  ON  SPECIFIC  COLUMNS  OF  B  " 

017452.  017452. 1.  1 

"CSCALE(NR) . THE  SCALINC  FACTORS  ON  SPECIFIC  COLUMNS  6f  C  " 

1.  1.  1.  1.  1.  1 

"NSCALE  (NC) . THE  SCALINC  FACTORS  ON  SPECIFIC  COLUMNS  OF  N  " 

1.1 

"K I SCALE  (NC) . THE  SCALINC  FACTORS  ON  SPECIFIC  COLUWS  OF  KI  " 

1.1 

"YSCALE  (NM)  - . THE  SCALING  ON  EACH  MEASUREMENT  " 

500.500.32,8 

"CNTRLSCALE  (NC)  -  -  -  THE  COMPUTER  SCALINC  ON  EACH  CONTROL." 

1024. 1024 

" YDSCALE  (NC) . THE  COMPUTER  SCALINC  ON  EACH  DESIRED  OUTPUT  " 

8.  32 

"YMAX(NM)  - .  THE  VECTOR  OF  MAX  VALUES  OF  THE  MEASUREMENTS  " 

20.  20.  20.  30 

"YDMAX  (NC)  . . THE  VECTOR  OF  MAX  VALUES  OF  THE  DESIRED  OUTPUTS 

30.  20 


SCALEM2  Computer  Program  Example  Data. 


"NR.NM.NC.TP" 

5  a. 3.  05 

“THE  MINIMAL  DISCRETE  COMPENSATOR  DYNAMICS  MATRIX.  PHI.  IS:” 

-O  78263.1  74302.0  90592.0  96006.0.72946 

"THE  MINIMAL  DISCRETE  COMPENSATOR  INPUT  MATRIX,  GAFMAMIN.  IS:" 

-O  00068.0  00024. -0  00074. -O  01160. 0  04977.0  01823,0  10458. -1  910236 
•O  00076. -0  00083.  -0  00088. -0  01252 , 0 . 04972 . 0  00800. 0  08074. -  1  750846 
0  00214, -O  00052.0  00084. -0.00071. -O  10784. -0.05930, -0  00479.0  370566 
-O  00048.0  00001.0  00156,0  00007,0  00701.0  00129.0  00276.0  007696 
-O  00356.0  00115,0.00006.0  00581, -0.67143.0  13067. -2.77471. -1  64643 
"THE  MINIMAL  COMPENSATOR  OUTPUT  MATRIX.  HMIN.  IS:" 

-1  79786. 1  82S97  1  00000. -O  16434. 1  000006 
O  10394.-0  03719.0  36894. 1  OOOOO, 0. 111846 
0  OOOOO  1  00000.0  79317.0.03578.0  09061 
"N  MATH  I X  (NC  X  NC)  COFWAND  DISTRIBUTION  MATRIX" 

•4  829E-02.-9  785E-03.-2  378E-036 
7  004E-03.-7  883E-02.2  098E-036 
-1  17SE-02.-3  420E-03. 5  886E-02 
"THE  KX  MATRIX" 

-  26522. -  19939. -1  2212 
"THE  KXI  MATRIX" 

-  0000864.-  0001028.-  17365 
"THE  KXD  MATRIX" 

-  094118. -  72424. -3  0 


"BSCALE  (NM)  - . THE  SCALING  FACTORS  ON  SPECIFIC  COLUMNS  OF  B  ” 

1.1  1.  017452.  017452.  017452.  017452.  017452 

"CSCALE(NR) -  THE  SCALINC  FACTORS  ON  SPECIFIC  COLUMNS  OF  C  " 

1  1.1  1.1 

"NSCALE(NC) -  THE  SCALING  FACTORS  ON  SPECIFIC  COLUMNS  OF  N  " 

1.1.1 

"YSCALE  (NM) .  THE  SCALINC  ON  EACH  MEASUREMENT." 

32.  32.  32. 500.  500.  500.  500.  500 

"CNTRLSCALE  (NC)  -  -  -  THE  COMPUTER  SCALING  ON  EACH  CONTROL." 

1024. 1024. 1024 

" YDSCALE (NC) . THE  COMPUTER  SCALINC  ON  EACH  DESIRED  OUTPUT." 

32. 32.32 

'YMAX(NM) . ---  THE  VECTOR  OF  MAX  VALUES  OF  THE  MEASUREMENTS." 

50.  50.  30.  10.  10.  30.  20.  30 

'YDMAX(NC) .  THE  VECTOR  OF  MAX  VALUES  OF  THE  DESIRED  OUTPUTS 

20. 20. 20 


FILMED 
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